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SUMMARY 


The experimental database established by this investigation of the flow in a large 
rectangular turning duct is of benchmark quality. The experimental Reynolds num- 
bers, Deans numbers and boundary layer characteristics are significantly different from 
previous benchmark curved-duct experimental parameters. The turbulent entrance 
boundary layer case, at R e( i = 328,000, provides an incompressible flowfield which ap- 
proaches real turbine blade cascade flow characteristics. Although incompressible, the 
turbulent boundary layer thickness of 5 to 10 percent of duct half-width provides a 
very challenging flowfield for CFD code development and validation. The experimental 
results present all three components of mean velocity and turbulence intensity for six 
measurement planes in the 90° turning duct. The experimental data is available on 
magnetic tape at the NASA-Lewis Research Center. 
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1. INTRODUCTION 


The advancement of turbine engine durability and efficiency depends upon design 
improvements in stage loading, operating temperatures and material properties. Im- 
proved fluid mechanics and heat transfer codes are key to design advancement of blade 
loading and temperature distribution. The curved duct flowfield results from the com- 
bined effects of compressibility, viscosity and heat transfer, and it is one of the most 
difficult three-dimensional flows to model accurately. A requirement for developing im- 
proved fluid dynamics codes is the benchmark experimental data base, which provides 
a calibration standard for development and validation of fluid mechanics codes. 

The investigation reported herein is the study of internal flow in a large rectan- 
gular turning duct, with a 25.4 cm (10 in) square cross section. The objective of this 
investigation is to provide a complete set of benchmark data for the flow within a large 
rectangular turning duct. These data are to be used to evaluate, verify, and ''or modify 
three-dimensional internal viscous flow models and computational codes for predicting 
the flow in turbine engines. The experimental and analytical efforts were conducted 
under a coordinated multitask contract 

The data acquired in the three-dimensional viscous flowfield includes pressure 
and velocity distributions for comparison with the flowfield computed using three- 
dimensional viscous flow computer codes. The measurements constituted a detailed 
study of the sidewall, endwall, and convex-sidewall/flat endwall corner regions of the 
duct, for isothermal, adiabatic flow. All measurements were taken at two Reynolds 
numbers that provided both laminar and fully turbulent boundary layer flows approach- 
ing the duct turn. The flow velocities for both Reynolds numbers were low enough to 
remain well within the incompressible Mach number range. 
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2. EXPERIMENTAL TASK 


Flow in curved ducts is representative of a wide range of real flowfields such as are 
found in turbomachinery. The evaluation of the complex three-dimensional flowfield in 
a turbine stage is complicated by the viscous effects which dominate the heat transfer 
distributions. A clear understanding of the three-dimensional boundary layer develop- 
ment in the curved duct is essential to the formulation and validation of computational 
fluid dynamics and heat transfer codes. This investigation expands the curved duct 
experimental work of Taylor, Whitelaw and Yianneskis (ref. 1) to higher Reynolds 
numbers and Dean numbers. The resulting experimental data base provides a CFD 
baseline data set that is closer to real turbine blade passage flowfields. In this section 
the experimental facility and instrumentation systems are described in detail. The 
design, fabrication and calibration of all experimental equipment was conducted with 
primary emphasis on data quality to assure benchmark results. 

2.1 Duct Flow Facility Description 

The experimental facility design features modular tunnel components which allow 
flow measurements every 15° in the 90° bend and at one duct width upstream and 
downstream of the bend. The 25.4 cm (10 in) square cross section tunnel is designed 
with a 13 to 1 area ratio bell mouth, contoured to provide uniform flow, and is pow- 
ered by an induced-draft, variable speed, variable pitch, six-bladed fan. The tunnel 
is designed for incompressible flow and will produce, based on variable inlet length. 
Reynolds numbers of 50 to 350 x 10 3 at the entrance of the 90° bend for tunnel veloc- 
ities of 4.0 to 24 m/sec (13 to 79 ft/sec). These two flow conditions represent laminar 
and fully turbulent boundary layer profiles at the entrance to the 90° bend. Figures 1 
and 2 show the details of the assembled duct. 

Inlet Transition The inlet transition section is designed to provide a uniform flow 
profile to the test section. The walls are made from formed metal with a beilmouth 
shape. The beilmouth has an overall length of 101.6 cm (40 inches), with a 91.44 cm 


- 4 - 



ORIGINAL PAGE IS 
OF. POOR QUALITY 



- 5 - 


Figure 1. Curved Duct Entrance Region and LV System 
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Figure 2. Assembled Curved Duct Test Section 



(36 inch) square entrance and a 25.4 cm (10 inch) square outlet for a contraction ratio, 
AJ Ao, of 13. The honeycomb panel acts as a flow straightener and smoother to provide 
a smooth, laminar core flow into the duct. An 8x8 mesh screen is stretched across the 
bellmouth inlet upstream and downstream of the honeycomb panels. The honeycomb 
and screens at the entrance, as well as the suction type air supply design are intended 
to provide the smoothest, most uniform core flow possible at the duct entrance. 

The bellmouth contour was selected on the basis of fluid dynamic calculations 
to provide an essentially uniform one-dimensional core flow with developing boundary 
layers entering the test section. Smooth flow acceleration is achieved by contouring 
the inlet to avoid large negative axial pressure gradients. An approximate calculation 
predicted the boundary layer thickness at the bellmouth exit of 1.25 mm to 2.5 mm 
(0.05 to 0.1 in). Therefore, the boundary layer at the start of the bend is mainly 
produced by flow over the straight duct extension sections, for both the laminar the 
turbulent boundary layer flow conditions 

The inlet extension section was designed to provide additional length upstream of 
the test section to allow development of a fully turbulent boundarv layer at the start 
of the bend This section is straight. 50. S cm (20 inches) long and has a 25 4 cm (10 
inch) square cross section Flanges are provided for ease of installation and removal, 
thus simplifying duct alterations to facilitate changing the flow characteristics. One 
of these interchangable sections is configured as an instrument segment, containing an 
optical access window module, an instrument probe module, and a 24.1 cm (9.5 inch) 
long extension module. This modular section can be interchanged or rearranged to 
permit inlet conditions to be measured at 25.4, 50.8. 76.2 or 101 6 cm (10. 20. 30 or 40 
inches) upstream of the bend entrance and exit conditions to be measured at 25 4 cm 
(10 inches) downstream of the bend exit. 

The instrument module includes the Laser Velocimeter access window as a separate 
piece. Construction of the instrument module is similar to the extension section, i.e. 
flanged and made from 6.35 mm (0.25 inch) thick aluminum plates. The basic difference 
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is the inclusion of the optical access window module and a traversing probe mounting 
panel which allows pressure measurements, temperature measurements, or hot-wire 
measurements, in the plane of the view window. Interchanging the instrument module 
with different panel modules permits flow conditions and profiles to be measured at 
the bend entrance, or at any of 5 evenly spaced stations in the bend. The view window 
module is constructed for application with either the curved instrument module or the 
duct-extension instrument segment. The window section was assembled from four (4), 
beveled 0.635 cm (0.25 inch) thick by 1.27 cm (0.5 inch) wide by 26.67 cm (10.25 inches) 
long flat glass panels mounted in an aluminum frame that bolts to the upstream end 
of the instrument module. The glass windowed module is arranged so that there is a 
clear view from corner-to-corner on all four walls of the test duct. Flat window width 
has been minimized within the operating constraints of the Laser Velocimeter system. 

Although flat windows were used along the convex and concave walls the resulting 
variation in wall curvature is insignificant. These 1.27 cm (0 5 in) wide flat panels 
were, by design installed so the window edges set on the arc of the concave and convex 
sidewalls, resulting at window center, in only a 0.004 cm (0.002 in) maximum variation 
from the desired uniform 45.7 cm (18 in) radius of the convex sidewall, and 0.003 cm 
(0.001 in) maximum variation from the desired 71 1 cm (28 in) uniform radius of the 
concave sidewall 

Induced Draft Fan, Exhaust Transition The exhaust transition section is designed to 
diffuse the flow and provide a uniform suction for moving the air through the test 
section. The exhaust section is 335.28 cm (132 inches) long and transitions from a 25.4 
cm (10 inches) square cross section at the downstream extension exit to an 82.55 cm 
(32.5 inch) diameter circular cross section designed to match the diameter of the fan 
mounted at the exit of this section. 

The selected, induced draft fan is a 6 bladed. adjustable pitch. 81.30 cm (32 in) 
diameter air propeller. The blades and hub are separate aluminum alloy castings, 
machine matched and balanced. The propeller is belt driven by a variable-speed. 
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electric motor mounted below the fan. Both the motor and propeller mounts are 
isolated from the duct supports, to inhibit any transfer of vibrations from the fan. 
Independent mounting arrangements were used to isolate exhaust transition section 
and the test section from fan induced vibrations and pressure fluctuations. Vibration 
isolation was provided at the flange connection between the exhaust transition section 
and downstream extension section, as well as between the exhaust transition section 
and the motor housing. 

Following facility construction, an extensive shakedown and calibration program 
was conducted to confirm flow quality in the inlet section and assure no flow separation 
in the diffuser section. Both laminar and turbulent inlet boundary layers were validated 
by hot-wire measurements. Flow uniformity in the inlet section was confirmed by 
pressure and laser velocimetry surveys. 

2.2 Instrumentation Description and Calibration 

The primary instrumentation was designed for non-intrusive flow measurements 
utilizing a three-component, laser velocimeter (LV) and wall static pressure gauges The 
LV utilizes two color beams and Bragg diffraction beam splitting/frequency shifting to 
separate the three simultaneous, orthogonal, velocity vector components. The LV signal 
processors determine the digital values of velocity from the seed particles crossing the 
laser beam probe volume. To improve and speed up digital data acquisition, the LV 
processors were designed around an S-100 bus Z-80 microprocessor which provides on- 
line, near-real time data reduction This on-line data was acquired and recorded for 
off-line detailed analysis. To qualify the measurements as "benchmark” data, the LV 
data was compared with both pitot probe and hot-wire anemometer measurements for 
flow conditions which permit comparisons. 

Total and static pressure distributions were measured at the selected measurement 
stations and along the duct walls. Total and static pressures measured at the exit of the 
duct inlet section were used to control the duct flow. Static pressure ports distributed 
along the centerlines of all four duct walls provided wall static pressure profiles along 
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the side and end walls of the duct as shown in Figure 3. 

Nominal bulk flow velocities of 6.0 and 20 m/sec. (20 and 66 fps) were selected 
to provide the laminar and turbulent boundary layer flows at the entrance to the turn. 
For these velocities, the pressure differentials to be measured are 0.10 and 2.5 cm of 
H 2 O (0.04 and 0.95 in. of H 2 0). To provide accurate measurements within this range, a 
Validyne DP103 differential pressure transducer for ultra low range was selected. This 
transducer has a combined accuracy of ± 0.5 percent of the calibrated span. The span 
can be altered by inserting different diaphragms. 

The full scale reading corresponding to the low pressure range measurements is 
1.0 cm in H 2 0, and the scale reading corresponding to the high pressure range was 
selected to be 4.0 cm // 2 O. Since a large number of static pressure measurements were 
necessary, two 48 port Scanivalve pressure scanners were used in conjunction with 
the Validyne transducers, to multiplex the static pressure readings. Therefore, all the 
pressure measurements were performed with three (3) similar pressure transducers, one 
for each scanivalve and one monitoring the inlet total to static pressure difference. 

Hot-wire measurements were made periodically in the flow for comparison with the 
LV velocity determinations. The hot-wires were calibrated in a steady, low turbulence 
flow' and the direct non-linearized output of the anemometer was recorded. From 
these data, statistically defined velocity measurements were determined and compared 
with LV measurements from the same flow, to ascertain possible biases in the LV 
measurements 

Laser Velocimeter System 

The laser velocimeter (LV) developed at UTSI and used on the curved duct flow- 
experiments may be described as two primary systems: (1) optical, and (2) electronic 
(ref 2.). The optical system consists of the laser, parallel beam color separator, two 
dimensional Bragg cell, lens train, lens train housing, photomultiplier tubes, and mill 
bed scanning system. The electronics system consists of Bragg cell drivers and power 
amplifiers, signal separator and mixer, simultaneous signal detector, signal measuring 
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Figure 3. Duct Static Pressure Port Location 
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electronics, data acquisition computer, and a scan position control system. The LV 
system was designed to measure three orthogonal components of velocity and to scan 
over a 25.4 x 25.4 cm cross section with a spatial resolution of approximately ± .25 
mm. 

A sketch of the optical system is shown in figure 4. A Lexel Corp. Model 95-2 
Ar 4 " laser was used in the LV system. This device produces approximately 2 watts of 
power in all line operation. It requires water cooling and 208 V, 3-phase, 30 amps/leg 
electrical service. The system used the 488 nm and 514.5 nm lines in LV operation. 
Guaranteed optical power for the 488 nm line is 0.7 watts and that for the 514.5 nm 
line is 0.8 watts. 

The color separator is a prism-mirror assembly designed at UTSI. Its function is 
to separate the input multicolor laser beam into a parallel set of beams which can be 
adjusted in terms of spatial separation and position relative to the input beam. The 
color separator consists of two glass dispersion prisms mounted such that the beams 
are incident on the prism surfaces at Brewster’s angle to minimize reflections. The first 
prism provides the initial color dispersion. The refracted dispersed beams are reflected 
by a first surface mirror to a second prism which further spreads the beams and reorients 
each beam such that it is parallel to the original input beam The intermediate mirror 
and prisms are adjustable in order to vary beam separation and output position relative 
to the initial beam direction. 

The parallel beams from the color separator are directed to a single two-dimensional 
Bragg cell which is used to split and frequency shift the parallel input beams. Each 
input beam is split into four beams of equal intensity originating from a common point. 
For a single Bragg cell each set of split beams appears to originate and propagate in a 
common reference plane For a two-dimensional Bragg cell, each set of split beams ap- 
pears to originate at the apex of a pyramid. The two-dimensional Bragg cell is a small 
water tank containing two x-cut quartz crystal transducers mounted at right angles 
with respect to each other. Each crystal is made to oscillate in an odd overtone 
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by a power amplifier referenced to a phase locked loop (PLL) oscillator. The frequency 
of each PLL is tunable in order that the transducer can be operated at its optimum 
frequency in a particular overtone. The output power of each power amplifier is ad- 
justable from 0-10 watts and allows control of relative beam intensities. The quartz 
transducers are cut for 5 MHz fundamental frequencies and are operated in third and 
fifth overtones to produce frequency shifts in the beams of 15 and 45 MHz. Frequency 
differences in orthogonal beam pairs are 15 and 45 MHz. Non-orthogonal beam pairs 
have frequency differences and sums of 30 and 60 MHz. Signal currents generated at 
these frequencies are removed by electronic filters after photodetection. 

After the beams are split, mirrors and beam stops are used to separate colors such 
that two orthogonal velocity components are measured with the 514.5 nm (green) beam 
and one component is measured with the 488 nm (blue) beam. All other split beams 
and colors are blocked. The green beams have velocity component identification and 
directions determined using the 15 and 45 MHz Bragg cell frequencies The blue beam 
utilizes a 15 MHz Bragg cell frequency. After the colors are separated, a dove prism is 
used to rotate the blue beams such that the velocity component, detected by the pair 
will be orthogonal to the component measured by the green beams. 

A lens was used to make the green and blue beams focus to a common volume 
of incidence. Transmission paths were matched to ensure that the probe volumes are 
approximately the same dimension. While the geometry of the volume is that of two 
intersecting cylinders, the volume will be contained within a sphere of approximately a 
200 jiim radius. Folding mirrors are used after the transmitting lens to assure sufficient 
adjustment for probe volume alignment. 

The entire optical assembly was mounted in an I-beam channel assembly for me- 
chanical stability. One side of the assembly was designed for quick removal in order 
to swing the system off the duct for alignment checks, etc. The channel assembly w'as 
mounted on a mill bed which had been modified to hold the channel and to producp 
a planar scan. The mill bed has positional resolution of 25 ^m and can be adjusted 
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manually or automatically with stepping motors under computer control. 

Positioning of the mill bed system at each of the measurement locations in the 
duct was accomplished by rolling the system on casters to the desired position. The 
mill bed was locked in position with stabilizing jack screws which raised the mill bed 
off the caster wheels and leveled the entire system. Zero coordinate reference for the 
computerized system was defined at the center position of the window assembly which 
was used at all stations. This position was determined by using the optical system with 
the blue and green beams to locate the window reference. 

LV Data Acquisition 

The electronics hardware associated with signal measurement and data acquisition 
is described in this section. Velocity was was acquired by three laser velocimeter burst 
counters. The fringe counts can be programmed from 2 to 255 which allows the operator 
to match the fringe count to the number of fringes in the LV probe volume Each burst 
counter is actually a combination of two time period counters, each of which performs an 
average signal period measurement. These period counters are started simultaneously, 
but one counter, called the short counter, stops after fewer fringe crossing than the 
second counter, called the long counter. That is. the short counter measures the signal 
time period averaged over .V s signal periods, and the long counter measures the signal 
time period averaged over N l signal periods. The significance of long and short is 
merely that Si > Ns The signal time period, ti. measured over Si signal periods 
and rs, measured over N$ signal periods are given by 


n-L ns 

TL = -Vi / a = Ssf..k 

in which / c jt is the clock frequency (250 MHz), ni is the number of clock pulses counted 
m time Sit l . and n$ is the number of clock pulses in time Ns's- 

The two measurements of the signal time period, ~i and t$ should, of course, agree 
because the measurements are performed during overlapping times on the same signal 
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The “aperiodicity” is defined by 


T L - TS 

e = 

tl 

Hard wired logic performs the aperiodicity test for each measurement. The operator 
can select a maximum aperiodicity of 1/128, 1/64, 1/32, or 1/16. Alternatively, the 
clock counts n i and n$ can be inputted to a computer which performs the aperiod- 
icity test in software, in which case the maximum aperiodicity limit is continuously 
adjustable. Experience has shown that tl, and t$ will either agree closely or disagree 
grossly. As a result, mean velocity and turbulence intensity measurements are only very 
weakly influenced by the precise value of the aperiodicity limit. Because the hardware 
aperiodicity test is much faster than the software test, one normally uses the former. 

Simultaneous signal detection in blue/green probe volume is required for turbulent 
data analysis. The signal time period measurement process is performed on each of 
the three velocity channels. Each burst counter outputs a pulse whose duration is 
.\ ' £ rt, called the ACTIVE pulse. The ACTIVE pulses from the three burst counters 
are passed to the simultaneity checker. Also, the aperiodicity fail pulse. AFAIL. from 
each burst counter is sent to the aperiodicity checker Only those data for which the 
three ACTIVE pulses were coincident for at least one instant in time, and for which an 
AFAIL pulse was not produced by any of the burst counters are passed to the computer 
which records the velocity data. 

The LV optical system was configured as a Bragg-shifted, two-color, velocimeter 
In essence, it is the orthogonal combination of two, two-component systems. Two 
PMT detectors were used and, hence, two UTSI LV signal receivers were required for 
frequency multiplexing two components of velocity from a single PMT signal. 

The signal splitter divides the total signal amplitude into two equal parts, each of 
which is passed to a signal processing channel consisting of a double balanced mixer. 
PLL local oscillator, low pass filter and burst counter LV signal processor. The signal 
channel responsible for processing signals produced by the 15.0 MHz Bragg transducer 
is called the 15.0 MHz channel, and the other signal channel is called the 45.0 MHz 
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channel. The LO (local oscillator) port of each mixer is driven by a PLL whose fre- 
quency is adjustable. For each signal frequency which appears in the mixer’s input (its 
radio frequency, or RF, port) the mixer’s output (its intermediate frequency, or IF, 
port) contains the sum of this frequency and the local oscillator frequency. However, 
the mixer output is passed through a low pass filter which lets through only those 
difference frequencies which fall within the passband of the filter. 

The signal frequency produced by a particle at rest, f 0 , is the difference between 
the frequency of the local oscillator, fio, and the Bragg excitation frequency, /#, 

fo = fto ~ /b 

A component of velocity for a moving particle is found from the difference between 
f 0 and the frequency of the output of the low pass filter, /. The latter is measured by 
a burst counter LV signal processor. (Actually, the burst processor measures r. the 
signal period averaged over some integer number of periods). Velocity was computed 
from the equation 

t’ = S(f 0 - /) 

in which 6 is the fringe spacing (i e . distance from one interference fringe to the next) 
The magnitude of a component of velocity was determined from the absolute value of 
fo - f. and the sign of the component of velocity was determined from the sign of -/. 

The LV translation table was controlled by the same computer which controlled 
the burst counters. A disk file of xyz positions was created for each desired flow profile 
and the data acquisition was automatic. In order to minimize operator error, the 
acquisition was automated as much as possible. 

LV Calibration 

Instrumentation and facility calibration was accomplished before the detailed flow- 
field measurements were obtained to assure benchmark quality data. The following 
paragraphs describe briefly the calibration standards and techniques for all instrumen- 
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tation. The calibrated instrumentation was first used to calibrate the facility (i.e. bulk 
velocity) and confirm core flow uniformity and boundary layer characteristics. 

The precision Validyne pressure transducers were calibrated before and after each 
data run with a precision slant manometer with full scale range of 13 mm H 2 0 and an 
accuracy of 0.1 mm H 2 0. These pressure transducers have a very linear response with 
demonstrated repeatability and accuracy of ±0.5 percent of calibrated span. 

Two techniques were used to calibrate the Laser velocimetry (LV) system for mean 
flow velocity. Then a single hot-wire system was used to check the mean turbulence 
intensity. A direct and absolute end to end calibration of the LV system was accom- 
plished using a spinning wheel device. A small target wire attached to the spinning 
wheel was rotated through each set of orthogonal fringes at a known velocity. Errors in 
knowledge of wheel radius and rotational speed were determined to be <0.5 percent of 
the 20 m/s calibration speed or 0.1 m/s. Tunnel velocity was also determined by pitot- 
static pressure measurement which was then compared with LV measurement. Due to 
the small pressure differentials associated with 6 and 20 m/s flows, the pitot-static de- 
termined velocity was at best a —1.0 percent measurement. During the experimental 
program, differences between pitot-static velocity and LV velocity measurements, at 
duct centerline, greater than ±1.5 percent were cause for recalibration 

Since the hot-wire anemometer is a well established measurement technique for 
turbulent flow properties measurement, a single hot-wire probe system was used to 
validate the LV measurements. Figures 5 and 6 show a comparison of mean velocity and 
turbulence intensity measurements taken by hot-wire and LV. In regions of the flowfield 
where high velocity gradients occur, the LV system will interpret the gradient produced 
velocity differences across the measurement volume as turbulence. This effect will 
result in erroneous measurements of pseudo-turbulence in thin laminar boundary layers. 
In the turbulent boundary layers and m the main duct flowfield. the LV measured 
turbulence intensity levels were 0.5 to 1.5 percent higher than the hot-wire levels. The 
significance of this difference will be discussed in section 4. 
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Figure 6. Hot Wire/Laser Velocimeter Comparison (Turbulent) 



LV Seeding 


The development and operation of a satisfactory flow seeding system was a difficult 
task which required several iterations. The particles must be small enough to track 
the flow accurately, but large enough to scatter sufficient light for accurate signal 
processing. The duct flow rates for this investigation were about one cubic meter 
per second which required a large amount of seed for an eight hour run. Due to the 
health hazards associated with most suitable solid particles such as aluminum oxide 
and titanium dioxide, they were rejected. 

Phenolic micro-balloons, which are used to reduce density in fiberglass, were found 
to be a good seed material. They have an equivalent particle size of 2 microns with 
good scattering properties. The micro-balloons were successfully seeded in both a 
water-alcohol slurry and in a fluidized bed air seeder. However the phenolic spheres 
were easily charged by static electricity and would collect on flow screens and duct 
walls. The operational problems of using the micro-balloons were greater than their 
advantages. 

Water droplets were selected for the comptete test program because of their good 
scattering characteristics, clean properties and low cost. Water droplet seed must be 
in the 2-5 micron size range to minimize particle lag problems. Water droplets as large 
as 10 microns followed the steady component of velocity with only a one percent error 
at mean velocities of 20 m/s due to the small velocity gradients in the curved duct. 
However accurate measurement of the low turbulence intensities associated with this 
project requires droplets of 2 microns or less. A water fog sprayer was used to generate 
particles in a large room from which the curved duct air was drawn. The larger particles 
were allowed to settle to the floow before reaching the inlet Thus a nominal 2 micron 
water droplet entered the duct. This technique worked well in the moderate relative 
humidity range 60 to 80 percent. L nder high relative humidity conditions the particles 
would grow and under low humidity conditions they would evaporate. Thus testing 
was restricted by weather conditions to moderate relative humidity days. 
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2.3 Duct Flowfield Measurements and Test Conditions 


Following checkout and calibration of the curved duct facility and instrumentation 
system, detailed flowfield surveys were taken with the three-component laser velocime- 
ter at six duct stations (inlet, 0°, 30°, 60°, 90° and exit). The inlet and exit stations 
were located one duct width upstream and one duct width downstream of the 90° 
curved section. Six hundred spatial points were surveyed at each of the six measure- 
ment stations for two Reynolds numbers corresponding to duct bulk velocities of 6 and 
20 m/s. The 6 m/s bulk flow case produced a laminar boundary layer at the entrance 
measurement station one duct width upstream of the 0° station. However, due to the 
physical size of the large turning duct, laminar flow could not be maintained through 
the 90° turn. A hot wire check of the exit station boundary layer confirmed a fully 
turbulent flow at bulk flow velocity of 6 m/s. Thus the “Low Reynolds Number” data 
has a laminar entrance boundary layer profile with transition to turbulent flow in the 
turn. The 20 m/s bulk flow case produce a fully turbulent entrance boundary layer 
which of course remains turbulent throughout the 90° turn. Thus the “high Reynolds 
number” data has fully turbulent boundary layers. At each measurement point approx- 
imately 300 data samples were taken by the LV system and processed to yield the mean 
velocity and its three orthogonal components along with average turbulence intensity 
and its three components. The turbulent cross-correlation terms are not presented in 
this report caused by the lack of statistical confidence due to the small sample size. 

To provide additional information on the flowfield, static wall pressures were mea- 
sured along the centerline of each duct wall. The wall pressure data shows the significant 
difference in inner and outer wall pressure gradients. In addition to the duct wall static 
pressures, the total pressure distribution was surveyed two duct widths upstream of the 
0 ° station. The total pressure was uniform outside of the wall boundary layer which 
was approximately one cm thick. The LV survey two widths upstream of the 0° station 
also confirmed a flat core velocity profile with thin boundary layer. 

For comparison, the experimental data obtained during this investigation was ob- 
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tained at much higher Reynolds numbers and Dean numbers than the data from NASA- 
CR-3367, reference 1. (See Table I.) The resulting thin boundary layers yield signif- 
icantly different crossflow velocity fields which provide a much more difficult baseline 
case for the computational fluid dynamics codes. 

Table I 

UTSI Investigation NASA-CR-3367 


R ed : 98,000 -+ 328,000 
De : 45,800 -> 152,800 


R ed : 790 — 40,000 
De : 368 18,700 
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3. ANALYTIC TASK 


To aid in defining the experimental program and analysis of results, and to de- 
fine the capabilities of state-of-the-art computer codes, this investigation contained a 
theoretical fluid mechanics task. A computer code was selected which was capable of 
solving the Navier-Stokes (N-S) equations with turbulence models for three-dimensional 
internal flow. The selected computational fluid dynamics (CFD) code was adapted to 
the experimental geometry and flow conditions of the curved duct. Validation of the 
code was accomplished by comparison of results with published experimental results 
and by making self-consistency checks of the computer results. Following validation of 
the selected CFD code, grid sensitivity studies were conducted to select an adequate 
grid spacing for the laminar and turbulent boundary layer calculations. Finally com- 
putations for the experimental curved duct configuration were made for comparisons 
with experimental data. 

3.1 Selection of 3-D CFD Code 

Several N-S codes were considered initially. Based on proven applicability to inter- 
nal flows similar to those in this investigation, two were evaluated in detail. The Briley- 
McDonald "MINT'’ code and the P. D. Thomas internal-flow code were both considered 
acceptable with relative advantages and disadvantages The Briley- McDonald code was 
judged more difficult to use and had much less available documentation. Also, the code 
was not available to us during the initial portion of the contract period. However, the 
“MINT” code had been run on curved-duct geometries similar to the configuration used 
in this investigation. The Thomas code, by contrast, has extensive documentation and 
had been applied to duct flows, although with different geometries from the present one. 
They differ mainly in the treatment of boundary conditions. Both codes are based on 
the Beam- Warming implicit, alternating-direction convergence scheme and are similar. 

For these reasons the P. D. Thomas internal flow' code was selected. It is fully 
elliptic with all inertial terms retained. Both cross stream viscous terms are retained, 
and only the streamwise viscous term is neglected. The Thomas code uses generalized 
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coordinates and was readily adapted to the curved duct geometry. Both code- validation 
and grid-sensitivity studies were conducted and are reported in the following sections. 

3.2 Code Adaptation and Validation 

The P. D. Thomas internal flow code was obtained from the NASA Langley Re- 
search Center and modified for use on the UTSI VAX 11/780 computer as well as 
on the, CRAY 1-S computer at NASA-Lewis Research Center. Two test cases were 
successfully run to demonstrate code baseline operation: laminar, flat-plate flow, and 
three-dimensional flow through a nozzle. Results were successfully obtained for nozzle 
throat Mach numbers of 1.0 and 0.5. The results were similar to those for test cases 
reported by Thomas (ref. 4). 

To validate the capability of the P. D. Thomas code extensive comparisons were 
made with the numerical results published by Briley and McDonald in reference 3 
and with the experimental data published by Whitelaw et. al. in reference 1 for flow 
in a square cross section turning duct. The comparisons were made for laminar flow 
to avoid the uncertanties involved in turbulence modeling The validation involved an 
assessment of the adequacy of the N-S solutions with practical computational grids, and 
an evaluation of the discretization errors. The basic comparison case was run at Mach 
number 0 1 and Reynolds number of 790, based on duct whdth. Although 20 percent 
fewer mesh cells than Briely and McDonald reported w r ere used for the comparison 
case, good agreement was found as shown in figures 7 and 8. Note that the comparison 
case selected has a large inlet boundary layer thickness of approximately 25 percent 
duct width. Both the P. D. Thomas and Briley-McDonald codes produced solutions 
that differed similarly from the experimental results. It can be seen that although the 
Briley-McDonald results used 69 percent more grid points in the crossflow direction, the 
results are very similar The fact that the peak in the experimental velocity profiles is 
not adequately resolved implies that both solutions are affected by discretization error 
with the grid sizes used. Thus a study of the influence of grid resolution was conducted, 
along with, a study of time step size and of smoothing functions for improved accuracy 
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Figure 8. Axial Velocity Profile Comparison, Z = 0, 0 = 90° 



and convergence. 


3.3 Grid Refinement Studies 

The grid selected for the validation case described above had 26 x 28 points in the 
cross stream plane and 17 streamwise planes. To investigate the effects of numerical 
truncation error, runs with different size grids were made. The selected grids had 
21 x 21 and 32 x 32 points in the crossflow plane and 17 and 33 down stream planes 
respectively. Figures 9 and 10 show only small differences in the results, which indicates 
a small truncation error. Also, the mesh was locally refined by a factor of three on the 
21 x 21 grid near the wall where a peak in the crossflow velocity occurs. No significant 
velocity profile changes resulted. For the entrance region profiles and Reynolds numbers 
evaluated, the resulting solutions were not sensitive to grid size or distribution. The 
discrepancy between the numerical and experimental results remains unresolved. The 
results do suggest that state-of-the-art N-S codes may have problems in resolving the 
detailed features of any complex flow. 

3 4 Numerical Results for Laminar and Turbulent Boundary Lavers 

Following extensive code adaptation, validation and refinement, flow-field calcula- 
tions were made for two Reynolds numbers corresponding to laminar and turbulent 
entrance boundary layers. Since the calculations were completed before the results of 
the experimental task were obtained, the cases were run at the same Reynolds numbers 
as reference 1. The laminar case was run at R e d = 790. and the turbulent case was run 
at R'd = 40,000. Both cases were run at Mach numbers of approximately 0.1. For the 
laminar case an inlet boundary layer thickness of 25 percent was selected, again to be 
compatible with reference 1. For the turbulent case, the algebraic eddy viscosity model 
contained in the Thomas code was used, and the sensitivity to initial boundary layer 
thickness was investigated. Comparative runs were made with 13 mm (0.5 in), 25 mm 
(1.0 in) and 51 mm (2 in) boundary-layer profiles and identical grid spacing. The 13 
mm and 25 mm initial boundary layer solutions were very similar, but the 51 mm 


- 28 - 





29 




Grid Study, y = 0.0, 0 = 
Reported Data 

Briley-McDonal d MINT Code 



Figure 10. Grid Refinement Study, Radial Velocity 




boundary layer solution was considerably different. These results may be explained, 
qualitatively, by considering vorticity effects in the flowfield. In all three cases, the 
vorticity entering the duct turn was about the same. In the first two cases, the initial 
vorticity was confined to a region near the duct wall, while in the third the initial 
vorticity was distributed over 40-percent of the duct half-height. The similarity of the 
two “thin” boundary layer solutions may indicate an axial numerical diffusion caused 
by insufficient grid resolution for the large gradients in the boundary layer, or may be 
a realistic physical result where the downstream flow depends on only the total inlet 
vorticity as long as it is confined to a thin region. Indeed, Thomas has even suggested 
that the number of grid points in the boundary-layer region be reduced to one. A check 
on mass conservation was made for both the laminar and turbulent flow calculations, 
axial mass flux was integrated over each cross section plane, the resulting mass flow 
variation was less than 0.5 percent. 

Figures II and 12 present the calculated axial velocity profiles for the curved duct 
exit station (6 - 90°) for the laminar and turbulent cases, respectively Note that 
both sets of velocity profiles are very similar and do not show significant influence of 
Reynolds number or boundary-layer tubulence. Although figures II and 12 agree well 
with Whitelaw's results from reference 1 for thick boundary layer flow at low Reynolds 
number, the next section will show less satisfactory agreement with experimental data 

Figures 13. 14. and 15 present the calculated crossflow results for axial stations 
0 = 30 °. 0 = 60 '. 9 = 90°. The resulting circulating flow pattern, which contains 
high crossflow velocities in the boundary layer, is dominated by the thickness of the 
inlet viscous layer. Lack of agreement with the experimental data presented in the 
next section can be attributed to insufficient grid resolution and unrealistically thick 
entrance boundary layers for the numerical solution initial conditions. 
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Figure 11. Axial Velocity Profile, Laminar, 
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Figure 12. Axial Velocity Profile, Turbulent 
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Figure 13. Crossflow Velocity Field, 
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Figure 15. Crossflow Velocity Field. 



4. PRESENTATION OF RESULTS 


The experimental data was acquired at two flow conditions, one at a low Reynolds 
number (98,000) corresponding to a duct bulk flow velocity of 6 m/s, and one at 
a high Reynolds number (328,000) corresponding to a duct bulk flow velocity of 20 
m/s. The low Reynolds number flow provided a thin laminar boundary layer at the 
entrance to the curved duct (6 = 0°). However, the Reynolds number was too high to 
maintain a laminar boundary layer throughout the 90° turn. This transitional flow was 
unavoidable because of the large physical dimensions of the curved duct and minimum 
stable operating velocity of 4 m/s. The high-Reynolds-number flow provided a thin 
turbulent boundary layer at the entrance of the curved duct. 

To assure that there were no extensive separated-flow regions in the curved duct, 
detailed wall static pressure surveys were conducted at conditions corresponding to the 
selected test cases. Inlet and exit total pressure surveys were also taken to evaluate duct 
flow quality. Typical experimental data plots are presented in the following section, 
and detailed data plots are presented in Appendix A and Appendix B. A complete 
set of data from the experiments is contained on magnetic tape at the NASA-Lewis 
Research Center and is listed in Appendix C and Appendix D 

4 . 1 S ummary of Experimental Results 

Wall Static Pressure Distribution The wall static pressure distributions from the cen- 
terline of the suction wall, pressure wall and top wall were reduced to pressure coef- 
ficient form (C p ). Figure 16 presents the static pressure coefficients for the low and 
high Reynolds number cases. The smoothness of the C p profiles indicates the quality 
of the duct construction and sensitivity of the precision pressure transducers Many 
wall static pressure surveys were taken before and during the laser veiocimeter data 
acquisition and all results were consistent. One of the dominant duct flow characteris- 
tics is clearly shown in the pressure coefficient plots. The large axial pressure gradients 
on both the suction and pressure walls induce the large crossflows in the low energy 
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boundary layer. Note that, in non-dimensional C p form, the pressure gradients are 
nearly identical for the low and high Reynolds number cases. The only difference in 
pressure distribution between the two cases presented in Figure 16 is the larger C p drop 
along the duct for the low Reynolds number case. This difference is perhaps a result 
of a thicker boundary layer for the low Reynolds number flow. 

Axial Velocity Profiles The primary component of velocity in the curved duct flowfield 
is the axial or x-component. It is the boundary layer generated by the axial velocity 
that is driven by the curved-duct pressure gradient to form the complex crossflows in 
the y and z directions. Thus, it is important to first consider the development of the 
axial velocity in the curved duct. The mean velocities presented in this section have 
been non-dimensionalized by the inlet bulk velocity (6 m/s or 20 m/s). To correct 
all data points to the same bulk velocity, the velocity at the tunnel inlet centerline 
(y = 0, z = 0) was utilized as a reference. Data were obtained by computer-controlled 
scan files containing twenty to thirty spatial data points. A zero-zero reference data 
point was taken at the beginning and end of each scan file to check for tunnel stability. 
The velocity components in each scan file were then non-dimensionalized by the average 
axial velocity of the two zero-zero data points. An average of all the center-line axial 
velocities from all scan files for one measurement station was then non-dimensionalized 
by bulk velocity, and this non-dimensional centerline velocity was used to correct each 
non-dimensional velocity to the bulk velocity reference. This somewhat complicated 
procedure corrected the entire data set, for one measurement station, for small changes 
in tunnel operating conditions. For those applications where actual velocity is required 
simply multiply the low- Reynolds- number mean velocities (U, V. W) by 6 m/s and the 
high-Reynolds-number mean velocities by 20 m/s. 

Figure 17 and 18 present a comparison ot low and high Reynolds number axial 
velocity profiles at the Z = 0 plane for all six measurement stations. Note that at 
two duct widths upstream of the 9=0 entrance to the 90° bend the velocity profile is 
uniform with a boundary layer thickness of 5 to 10 percent of the duct half-width. The 
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b. High Reynolds Number Flow 

Figure 16. Wall Static Pressure Distribution 


- 39 - 








© 


ID 

O 


© 

© 


ID 

Q 

i 


© 

© 

© 

I! 

P-J 


X 

Cl 

ID 


m 

u 

c 

iTI 

Ml 

© 

r-H 

iTI 

•Q 

iTI 

LL 


©• 


iTI 


U .= U. 


i i_i 


O 

O 


II 


- 40 - 


Figure 17. Experimental Axial Velocity, Low R e d, Z 
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Figure 18 . Experimental Axial Velocity High R e d,Z = OiO 



velocity gradient between the suction wall to the pressure wall is well established at 
the 0 = 30° station, and remains nearly constant until suction-wall pressure gradients 
retard the flow between 0 — 60° and 0 = 90° . 

To illustrate the variation of axial velocity in the 2 -direction, Figure 19 compares 
the low and high-Reynolds number axial velocities for the 0 = 30° station. Six 2 -planes 
are shown from the duct center plane down to 4 mm from the bottom wall. Only the 
bottom scan shows any influence of the 2 -direction velocity gradient. This data has 
not been non-dimensionalized thus there is some scatter in the core velocity profiles. 
There is, of course, a much greater 2 -variation in the axial velocity at stations 0 = 60° 
and 0 = 90°. See Appendix A and Appendix B for additional axial velocity profiles. 

Crossflow Velocity Field The crossflow movement of the pressure gradient driven bound- 
ary layer has a very important influence on convective heat transfer. Thus the devel- 
opment of the crossflow velocity field in the curved duct was a major area of interest 
in this investigation. The crossflow vector plot provides the best presentation and eas- 
iest understanding of this flowfield. Figures 20 and 21 show the development of the 
low-Reynolds-number crossflow for the 0°,30 o ,60°. and 90 J measurement planes. At 
the 0 — 0 entrance to the 90° turn, the suction-wall and pressure-wall gradients have 
already influenced the crossflow, which results in a general flow toward the suction wall. 
The streamline curvature and radial acceleration effects have already reversed the gen- 
eral crossflow direction at 0 = 30° except for the top and bottom wall boundary layers, 
which continue to accelerate toward the suction wall. At the 0 — 60° plane the classic, 
curved-duct crossflow has been established. The magnitude of the crossflow velocities 
outside of the boundary layer are generally less than 5 percent of the axial velocity. 
However, in the developing boundary layer along the top, bottom and suction-walls 
crossflow, velocities of 20 percent magnitude are found. The crossflow velocities reach 
a maximum value at the 0 — 90° plane. The developing boundary layer has been moved 
by the strong crossflow and collects near the center of the suction wall in a region of 
high turbulence and low axial velocity. Downstream of the 90° plane, the 
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MEAN FLOW VELOCITY PROFILES 
LOW REYNOLDS NUMBER CASE - STATION 6 <0=30 DEG) 



Figure 19. Experimental Axial Velocity Profiles, 




STATION 4, 0° CROSSFLOW 


U k = 6 M/SEC 



U h = 6 M/SEC 



Figure 20. Experimental Crossflow Velocity, 6 = 0°, 6 = 30° 
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STATION 8, 60° CROSSFLOW 


U b = 6 M/SEC 



STATION 10, 90° CROSSFLOW 


U k = 6 M/SEC 



Figure 21. Experimental Crossflow Velocity, 0 = 60°, 6 = 90° 
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strong crossflow toward the pressure wall decays rapidly but the circulation pattern 
near the suction wall persists. 

The non-dimensional crossflow field for the high-Reynolds-number case is nearly 
identical to the low-Reynolds-number flow described above. 

Turbulence Intensity Distribution The mean turbulence intensity measurements pre- 
sented in this section were made with the 3-D laser velocimetry system. Average 
sample size per data point was 300 which is considered small for turbulence measure- 
ments. (This sample-size limit was imposed by the available microcomputer system.) 
A single-wire, hot-wire anemometer was used to check the LV turbulence intensity. 
In general the LV measurements of intensity were I percent higher than the hot-wire 
measurements. Thus the turbulence intensity values presented in this report should be 
considered high by l percent with an accuracy of xl percent. 

Figure 22 compares the low-Reynolds-number duct centerline turbulence intensity 
with the high-Reynolds-number values. The turbulence distribution is nearly identical 
for the two flow conditions. The core-flow turbulence intensity downstream of the inlet 
screens was measured at 2-4 percent which is a reasonable level The boundary layer 
intensity levels were measured at 6-10 percent which is also typical for thin turbulent 
boundary layers. Note that all turbulent intensities were non-dimensionalized with 
duct bulk velocity. 

4 2 Analys is and Evaluation of Results 

The experimental database established by this investigation is of benchmark qual- 
ity, and the experimental conditions are representative of actual turbine cascade flows. 
The experimental Reynolds numbers, Deans numbers and boundary-layer character- 
istics are significantly different from previous benchmark curved-duct experimental 
parameters Although incompressible, the turbulent flow Reynolds number of 328.000. 
based on duct width, combined with a 5 to 10 percent boundary layer thickness based 
on duct half-width, provides a very challenging flowfield for CFD code development 
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LOW REYNOLDS TURBULENT INTENSITY DISTRIBUTIONS 
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HIGH REYNOLDS TURBULENT INTENSITY DISTRIBUTIONS 
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Figure 22. Turbulence Intensity Comparison, Z — OiO 
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and validation. The experimental conditions produce large velocity gradients in the 
axial and crossflow directions in the relatively thin boundary layer regions. A CFD 
code which provides good results for low-Reynolds-number, thick-boundary-layer flows 
may not accurately solve the realistic flowfield benchmarked in this investigation. 

Since the numerical results were computed before the experimental results were 
obtained, actual test parameters and initial-boundary layer conditions were not uti- 
lized in the numerical investigation. Direct numerical-to-experimental comparison is 
not possible, however, the sensitivity study conducted with the CFD code highlights 
several areas of concern. Grid spacing refinement and improved turbulence modeling 
are the two recommended improvement areas for making the CFD code results agree 
with experimental results. The grid-refinement studies conducted in this investigation 
showed only a very small influence on the results, which were in poor agreement with 
experimental data. The addition of a simple turbulence model to the CFD computa- 
tions also showed small sensitivity in the results. Thus this investigation has produced 
an experimental database for a relatively simple flowfield which provides a significant 
challenge for CFD code development 
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5. CONCLUSIONS AND RECOMMENDATIONS 


This investigation complements the results of Taylor, Whitelaw and Yianneskis, 
ref. 1, and extends the experimental database to higher Reynolds number and thin- 
ner entrance boundary layers. The 5% to 10% thick boundary layers, based on duct 
half-width, resulted in a large region of near-potential flow in the duct core surrounded 
by developing boundary layers with large crossflows. The turbulent entrance bound- 
ary layer case at R e( i = 328,000 provides an incompressible flowfield which approaches 
real turbine blade cascade characteristics. The experimental results of this investiga- 
tion provide a challenging benchmark database for computational fluid dynamics code 
development. 

Based on the CFD studies conducted as a part of this investigation the following 
recommendations are offered. 

1. The grid spacing that was successful for calculating laminar flow with thick bound- 
ary layer will not be adequate for the thin boundary layer case. 

2. The thin turbulent boundary layer flowfield will require improvements in both grid 
resolution turbulence modeling, and perhaps in the computational algorithms 

3. Particular attention should be placed on the distribution of axial calculation planes 
in regions of high axial pressure gradient The axial velocity profile must be 
accurately calculated to assure proper development of the crossflow field 

4. The importance of the turbulence model and its influence on the CFD results 
was not determined in this investigation. However, the database established will 
allow the future investigation of the inadequacies of locally dependent turbulence 
models. 
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LIST OF SYMBOLS 


A Flow area 

C p Static pressure coefficient - 

Dh,d Duct hydraulic diameter (width, height) 
D e Deans number - (d/2r c )^ R e 4 

f Signal frequency 

N Fringe crossing count 

n Number of clock pulses 

P Pressure 

q Duct dynamic pressure (inlet) 

Red Reynolds number based on d 
R t Duct inner wall radius 

R 0 Duct outer wall radius 

r, : Mean radius of curvature 

U Non-dimensional velocity (x-direction) 

u x-velocitv component 

V Non-dimensional velocity (y-direction) 

v y-velocit\ component 

W Non-dimensional velocity (2-direction) 

w 2- velocity component 

x Axial coordinate (main flow direction) 

y Radial coordinate 

2 Vertical coordinate 

5 Fringe spacing 

e Aperiodicitv 

a Turbulence intensity 

r Fringe period 
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APPENDIX A 

Low Reynolds Number Data R e d = 98,000 


Figure A-l 
Figure A-2 
Figure A- 3 
Figure A- 4 
Figure A-5 
Figure A-6 
Figure A-7 
Figure A-8 
Figure A-9 
Figure A- 10 
Figure A-ll 
Figure A-12 


Wall Static Pressure Distribution 
Axial Velocity Profile, Z = 0.00 
Axial Velocity Profile, Z = —0.50 
Axial Velocity Profile, Z = —0.95 
Crossflow Velocity Field, 0 = 0° 
Crossflow Velocity Field, 6 = 30° 
Crossflow Velocity Field, 6 = 60° 
Crossflow Velocity Field, 0 = 90° 
Crossflow Velocity Field, Exit ■+ ID 
Turbulence Intensity, Z = 0.00 
Turbulence Intensity, Z = -0.30 
Turbulence Intensity, Z = —0.60 
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Figure A-l Wall Static Pressure Distribution 



LOW REYNOLDS MEAN FLOW UELOCITY PROFILES 


© © 
•ULDUL 3 • 
CM LU UJ UJ tU •»* 
I DDuD + 



© 

in 

© 

in 

© 

CM 

■H 


© 

© 


© X •*” AJ i-H 

O O £ CL 

3 \ O 



A- 2 


Figure A-2 Axial Velocity Profile, Z = 0.00 



Re< H > = 9 . 94E4 Del P=0.0S4" H20 

De< H >=4 . 63E4 IK B > =6 0 m-'s 


OX<J □ *<> 


© 

% 

© 

3 

'w' 

:> 

^•<r 
s \ 
£ A 
A 

Tf© 

in - 

OJ© 

V V 
:> 
C <£ 
v 

3 

© -W- 
/a 

ii in 
13 

ii 

QO 





LOW REYNOLDS MEAN FLOW UELOCITY PROFILES 


© © 
• UJ i J i J i JJ • 
CM UJ LU UJ LU 
I Cl Cl Cl Cl -f 



Q 

in 

© 

in 

© 

• 

• 

• 

• 

• 

CM 

-H 

^4 

© 

© 


^ X •*- iTJ •— i 

U 0 E CL 

D \ U 



A- 4 


Figure A-4 Axial Velocity Profile, Z = -0.95 
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Figure A-7 Crossflow Velocity Field. 6 = 60' 
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Figure A-9 Crossflow Velocity Field, E: 
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Figure A-10 Turbulence Iu tensity, Z — 0.00 
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Figure A-ll Turbulence Intensity, Z, ■= -0.30 
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Figure A-12 Turbulence Intensity, Z = -0.60 



APPENDIX B 

High Reynolds Number Data R e d = 328, 000 


Figure B-l 
Figure B-2 
Figure B-3 
Figure B-4 
Figure B-5 
Figure B-6 
Figure B-7 
Figure B-8 


Wall Static Pressure Distribution 
Axial Velocity Profile, Z = 0.00 
Axial Velocity Profile, Z = —0.50 
Axial Velocity Profile, Z = —0.95 
Turbulence Intensity, Z = 0.00 
Turbulence Intensity, Z = —0.30 
Turbulence Intensity, Z = —0.60 
Turbulence Intensity, Z — -0.80 
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Figure B-l Wall Static Pressure Distribution 
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Figure B-2 Axial Velocity Profile, Z = 0.00 
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Figure B-3 Axial Velocity Profile, Z = —0.50 
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Figure B-4 Axial Velocity Profile, Z = -0.95 
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Figure B-5 Turbulence Intensity, Z = 0.00 
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Figure B-6 Turbulence Intensity, Z = —0.30 
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Figure B-7 Turbulence Intensity, Z — —0.60 
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Figure B-8 Turbulence Intensity, Z = —0.80 


APPENDIX C 

Low Reynolds Number Tabulated Data 


Table C-l 
Table C-2 
Table C-3 
Table C-4 
Table C-5 
Table C-6 


Station 2, X(H) = 2.0 Entrance 
Station 4, 0 = 0° 

Station 6, 6 = 30° 

Station 8, 6 = 60° 

Station 10, 6 = 90° 

Station 11, X(H) = 1.0 Exit 


The tunnel fcoordinate positions have been non-dimensionalized on duct half-width 
Dh/z. All velocities have been non-dimensionalized on bulk velocity (6 m/s). 



Table C-l, Station 2, X(H) = - 2.0 Entrance 
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-0.084 

O. 044 

O. 045 

0 . 037 

0. 04u 

-0. 93 

0. .10 

\ . is 

1.25 

0. ft 02 

-ft. 072 

0. 045 

1.1 . 045 

0 . 035 

ft. ft 45 

-0. 93 

ft. 2ft 

1 . 34 

1 . 23 

O. urn 

-ft. US 9 

ft. ft .7 4 

(J. 034 

a . 030 

0. 03 7 

-0 . 93 

o . 1 V) 

1 . ?L, 

1 . 26 

Cl. ocill 

- 0.049 

ft . ft 35 

o. 0 35 

0. 039 

ft . ( >35 

“0. 73 

» i . Oft 

1 . 36 

1.26 

ft . oos 

-0.042 

O. 049 

0. 049 

ft . 034 

( > . u5 ? 

U . 9 3 

-1 1 . 1 Cl 

1.3'. 

1.7.1 

( ) . Of if ) 

- ft . r»:;s 

o. o34 

ft. ft 34 

O. 072 

C) . ( 1 • 5 

-0. 73 

— u. 2< ) 

1 . 3r> 

l . IS 

ft. oft *5 

-0.027 

ft . • r s 

0.035 

0. 075 

f > . 1 ) ". 7 

-0 . 9 3 

0. *V> 

1 . .5 

1 . 37. 

ft. ftftS 

-ft. ft 25 

0. ( >4S 

0 . 043 

0. 1)34 

P (>44 

- o. 73 

- n . 4o 

i . .3 

1 . _ ’. 

o . ft i 

- ft . ft l s 

(>. <>55 

ft. ftS5 

n.037 

0.(14 1 

-o. ? 

-i So 

1 . 3 1 

1 . 71 

i.i . 1 1 1 7 

ft . < /ft 7. 

o. on 5 

<1. Cm5 

Cl . 0 14 

(1.04 7 

- ft. V3 

-o. ,S» 1 

» . 1 / 

1.16 

n. 02 7 

U. ft _ 2 

0.001 

O . OB 2 

0.042 

1 1. 054 

-li. 93 

ft. 7ft 

1 . ( »B 

1 . 07 

0.03 7 

ft. (.‘34 

0.097 

II. 07/ 

O. 057 

0.1)57 

-i 1 . 7 

-( > . m > 

0. 9.1 

o. 7 7 

ft. i)7.5 

O. 045 

0 . 11 7 

0.11 9 

0. 070 

1 1 . 069 

9. 75 

o. it' i 

i .in 

1 . (|C| 

ft. Oft 3 

-ii. 177 

Cl. 07 > 

U . ft '/ 5 

0. 058 

o. u5 r 

-o. 75 

ft. ?o 

1 . Ol 1 

0. 99 

11.01 0 

-0.073 

o . 23 n 

ft . 253 

0.054 

0.1 >59 

- o. > ,u , 

0. 61 1 

l. In 

1 . 09 

-Ci. do? 

-ft. ft?5 

1 1 . OV 6 

n. i.i 78 

0. 049 

0. 056 

o. vs 

ft. SO 

1.17 

1.17 

( >. i >07 

o. 051 

ft. \>n ’ 

ft .1)113 

0.042 

ft. n 4B 

o. 75 

ft. 4o 

1 . 1 a 

1 . 18 

ft. OftO 

-O. 065 

ft . < >8ft 

O. 08' i 

U.II4 3 

" . I '4 6 

- '■> . 95 

• > . "’ft 

1.19 

1.19 

- l’l. 005 

-0.061 

o . i >n 

0 . 0(15 

0 . 039 

ft. ft" i 

-ft. V» 

ft . ' / 

1 . 1 

1.31 

-',i. mi/ 

* ft. ft56 

ft. UV- 

u. o9.' 

ft. 039 

ii. i '44 

~ ( > . ?S 

o. 1 o 

1 . ' 

1 . 33 

- ft . 005 

-ft. (>3o 

(./ . ( >fl 2 

ft. 1 Ul'. 

1 1. 034 

V. u 

-ft. 95 

( ‘ . ftt i 

1.31 

1.31 

( > . Oi >0 

-ft . 01 0 

ft. ft 75 

ft. 07 5 

0 . 036 

ft . < > *v 

- < i . 95 

--D. in 

1 . 33 

1 . 23 

—O. i if) 2 

- ft. u3? 

Ci. 070 

0. 1)70 

0.032 

( > . (> Z 4 

-ft. 95 

* Vo 

1 . „ 

i . 

0. OuS 

“U. 0 34 

Cl. (i£|3 

o . oar 

0. ft 36 

o 04 4 

•-o. VS 

-ft. r .o 

1 . ‘ 3 

1.22 

ft. fn >7* 

-0 . 020 

0. o/O 

Ci. 1 1 70 

0. 037 

ft. <>34 

I ) . 95 

“ft. 4ft 

i . ta 

i . i a 

0 . uft* ) 

-0.005 

( i . 01 1'i 

0 . OR3 

Cl. Ci 4 4 

I1.D41 

“O . vs 

--ft. S' > 

1.14 

1.13 

-t ft ft( » ? 

ft - ft 22 

0.087 

o . on 7 

O. 051 

0 . ft5 1 

-1 '. 75 

-Ci. 6i) 

1 . ' »5 

1 . t>5 

O. 007 

0. ' i4fa 

ft. 1 (>2 

O. M2 

0- V 5 V 

U. ft 39 

-o. 7*5 

-ft. 7o 

ft . 9 7 

11. 96 

O. 024 

ft . ( >6R 

O. 107 

0. 107 

Cl .065 

0 . 1 >66 

-U. 95 

~u. tin 

ii. 60 

o. 66 

0. Ol 7 

0.024 

o. 425 

0.45O 

0 ,u61 

( > . (>5 j 

-O . ?6 

ft. Oft 

o. 74 

0. V? 

-O. 014 

-O. 164 

0. U75 

Cl. 075 

Cl. 054 

ft. < >62 

-O. 9t> 

ft. /ft 

1 . ftO 

O. 99 

o. 007 

-O. 1 00 

0 . 1 09 

(i.lll 

1 1 . 066 

II. 1170 

-ii. 76 

. fc>ft 

1 . 05 

1 . 04 

- «>. Oo 7 

-0. 1 17 

0. to5 

O. 109 

0 . ft 55 

ft ftto? 

-0.96 

ft. S' > 

i.ii 

1 . 11 

-o. Oln 

-0. 098 

Cl. 104 

ii. 1 05 

o.o52 

ft. ft 5 7 

~U . V^J 

ft . 4 ft 

1.16 

t . 16 

-n. 007 

-0.078 

1 1 . 075 

0. 095 

0. 047 

ft. f>5< > 

- 0.96 

ft. .u 

1.17 

1 . Ill 

- ft . ft J 2 

-O.U71 

n. 098 

o. loO 

Cl. 04 8 

0. ft 4(3 

-0.^6 

ft. ; ft 

1 . 1 H 

i . m 

-o. 009 

-0. ( >55 

0 .(>92 

ft. 093 

0. 043 

U. 041 

-o. 96 

(.». i'» 

1 . 1 

1 . ?0 

-O. 0o7 

- 0.057 

o . i oi i 

O. 102 

O . ft 55 

0. i>4 n 

-ft. 7 6 

ft. Oft 

1.3’. 

1 . 23 

-0. 002 

-0. ft 2ft 

o. (>; i 

0. 07 1 

O . Cl 20 

ft. ft? ? 

- o . V 6 

“ 1 ft I o 

1 . 

1 . .2 

ft. ftftft 

-0. 028 

Cl. 069 

0 .0/1 

0 . 0 3 1 

(> . ft :.3 

-o. 76 

- ft . Ill 

1.31 

1 . 21 

--O . Oft 2 

-ft . 0?2 

U. 078 

0. i>79 

o. 035 

0.1)71 

- ft. 7<s 

-•ft ' t'» 

j . : i 

1.71 

-O.i >u3 

-( > . f>0 3 

ii. 0/8 

O. 1)73 

o. or. a 

0.1*2 

-ft. /<•» 

- ft . |o 

1.14 

1.14 

—»j . Oi >5 

ft.012 

0. 1 "2 

U. 1<>4 

0. 045 

U. ' 14 

-o. >/» 

-ft. S ( > 

1.11 

1 . 1C) 

o. 003 

V. 038 

0.1 >97 

0. 09 7 

0. 059 

o. 05» > 

-ii. 96 

-0 . 60 

1.01 

1 . U 1 

O. 007 

0 . 059 

O. 124 

o. 124 

0. 062 

0. 059 

-o. 96 

-t l . ?o 

( > . Vi ) 

Cl. B? 

0.014 

0. Ci 90 

O. 130 

O. 1 S3 

0 . 07 1 

li. 07? 

-0.96 

-0. Bo 

ii. 71 

0 . 90 

0 . 0 1 0 

ft. 026 

u. 1 12 

0.114 

Cl. 069 

0. 0/14 
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y 

z 

Urms 

U 

_ V 

w 

U rms 

u' 

v' 

VI 1 

O. So 

0. 00 

0 . 03 

0.83 

0. 003 

-0. 002 

0. 068 

0. 070 

0 . 06 1 

0 . 055 

0 . 70 

o. oo 

u.as 

0.87 

0. 0o2 

-0.005 

0. 060 

0.061 

0 . 055 

0. 042 

0 . AO 

0 . 00 

l>. 91 

0.90 

0.010 

-O. 007 

0.048 

0.050 

0.050 

17. 04? 

l.l. 50 

0. 00 

0. 93 

0.93 

0.010 

-0.010 

0.048 

0. 048 

0. 042 

0 . 030 

0.40 

0.00 

0.95 

0.94 

0 . ( >20 

-0.017 

O. 045 

0. 045 

0.O3B 

0. 038 

0. 30 

0. 00 

0.97 

0.97 

0. 028 

-0.012 

O. 043 

0.043 

0. 038 

0. 038 

0. 20 

0 . oo 

0.99 

O. 99 

0. 020 

-O. 01B 

0.045 

(>. 045 

O. 043 

0. 035 

i). 10 

0. OO 

1 . Ol 

l.til 

0. 038 

-O. 008 

0. 040 

O. 042 

0. 038 

0. 040 

o. 

0. 00 

1.03 

1 . 03 

0. 040 

-0. 01 ft 

0. 04 2 

O. 043 

0. 040 

0. 04 (7 

-0. 10 

* » . OO 

1 . 116 

1 . 06 

0 . 040 

-0. 012 

0. 042 

0. 042 

O. 040 

0 . 057 

-l). ?ll 

0. Ou 

1 . 07 

1.06 

0 . 050 

—0.013 

0 . 045 

o. 045 

O. 040 

ii. 050 

”0. 30 

0 . 0* » 

1.11 

1.11 

0 . OS l 

-0.015 

0.042 

o. 043 

O. 037 

0.037 

-ii. 4u 

ci. t >0 

1 . 12 

1. 1 l 

U . 06 1 

-0. 015 

0.043 

ft .04 3 

0.040 

0. 035 

-O. 5u 

o. oo 

1 . 15 

1. 14 

0.070 

-o. 017 

0.043 

0. 043 

0. 037 

0 . u3 7 

-<i.55 

o . u 1 1 

1.16 

1. 16 

0.071 

-0.020 

0 . 050 

ci. 05u 

0.042 

O. 040 

-0. 60 

0. Oo 

1.16 

1. 16 

0. 075 

-0.012 

0.045 

0. 047 

O. 043 

O. u40 

-O. 65 

0 . Oc » 

i . ta 

1.17 

0 . 070 

-0 . 008 

0. 043 

0.043 

0.035 

( i. 03 7 

-0. 70 

0. 00 

1 . 20 

1 . 19 

0.063 

-0.010 

0.048 

0.048 

0.037 

(i. 04 2 

-0.75 

o» oo 

1.71 

1.21 

0.076 

-0. 003 

0. 047 

O. 047 

0.037 

0.(7 38 

-0 . BO 

0.00 

1 . 22 

1.21 

0 . 065 

-0 . 002 

0 . 050 

0 . 05 1 

0. 040 

0. 038 

- 0 . BO 

-0. 1 O 

1.23 

1 . 22 

0.670 

ft .012 

0. 047 

0. i >47 

0. 038 

U . 031-1 

-0 . 75 

-o. to 

1 . 2u 

1 . 20 

O. 073 

-0. 007 

0. 043 

0. u50 

0. 035 

0 . 04 3 

-it. 70 

“0. to 

1.19 

1.19 

0.071 

-0. 013 

o. 047 

0 . 04 7 

(i. 037 

(1.1)43 

-0. 65 

— 0. 10 

1.19 

1. IB 

0.075 

-0. 015 

ft . 04 3 

0. 043 

0. 042 

0. 040 

-V. 61 1 

-o. to 

1.16 

1.16 

0.070 

- ii. 022 

0. 048 

o. 048 

0.042 

0.0 AO 

-0 . 55 

~U. to 

l . 15 

1.15 

O. 063 

- 0.018 

ft. 04 7 

0.(14 7 

0 . 038 

0.O37 

-O. 50 

- U. l V 

1.14 

1.13 

0. 068 

-ii. 028 

ll. i>4 7 

0. i>47 

0. (743 

U . 0 38 

-(I. 40 

-0. lo 

1.13 

1.12 

O. 056 

-O. 022 

O.045 

0. 045 

0.043 

1 1 . 035 

•• o . r.t > 

-0 . 1 o 

1 . 09 

1 . 09 

0.04 7 

-0.018 

0 . 043 

0.043 

0. 045 

O. 0 33 

-0. 20 

-0. to 

1 . OB 

1 . OB 

0. 047 

-0 .018 

0.043 

0. 045 

0. 042 

( i . 040 

-II. 10 

•ll. JO 

1 . 06 

1 . 06 

0. 045 

-(). 01 7 

o. 045 

o. 045 

O. 038 

0.053 

1 ) . 00 

- 0 . 1 c ) 

1.117. 

1 . 0.7 

0 . 035 

-IJ. 008 

(i .04 7 

0. 047 

0.042 

O. 0 33 

0 . 1 O 

- o . to 

1 .OJ 

1 .Ol 

U. 030 

-0. 0OB 

o. 040 

0.04? 

0.038 

1 1 . 035 

11.20 

-o. lo 

l . ou 

1 . 00 

O. 030 

-O.013 

0.037 

0. 038 

o. 047 

u, 1 1 3u 

II. 'll 

o. lo 

O. 97 

o. 9 7 

0. 027 

- 0 .013 

0. 042 

o. 04? 

O. 058 

o.o ’,5 

0. 40 

-0. J o 

1 1 . 95 

O. 94 

( >. 020 

-0. 01 2 

ii .040 

U. 040 

0.045 

(7. 035 

O. 50 

0. to 

O. 94 

0. 94 

0. 020 

- 'i . uon 

O. 040 

0. 040 

0. 043 

ii.o3. 

ii. 60 

-0. ID 

li. 92 

0. 92 

o. 005 

-0 . 008 

0 .07 3 

O. 033 

O. 055 

O. ii '.O 

o. 70 

C» . 1 o 

0.89 

0.1:17 

0. 0u2 

-O . OoQ 

i > . 055 

0. 055 

0.067 

0.045 

o. no 

"0. io 

o. ns 

0.05 

O. 003 

(i . 003 

0.067 

0 .06 7 

0.062 

0. <>43 

0.1 III 

o« J*u 

0. S3 

0.8.-; 

0.0 to 

O.ooS 

0. 065 

O. 065 

0. 06? 

ll. 04/ 

O. 70 

- o . .*o 

u.as 

O.HB 

0.013 

-0 . < 107 

0. 065 

0.067 

0.055 

0. 043 

0. nO 

“O . ‘U 

O. 91 

0.90 

(1.012 

-0.018 

0. 052 

0 . 052 

0.04 8 

0 . 0 - 

1 1 . 50 

~ (.» . r< > 

0.97 

0.9 2 

0. Ol 2 

-0.022 

O. 043 

0. 050 

0.047 

0.035 

0 , 40 

-o. : i o 

ll. 96 

0 . 96 

O. 032 

-0.015 

n. 045 

(>. 045 

0. 047 

O. 035 

O. 50 

—0. JO 

0.93 

0.97 

O . « >3 3 

-ft .017 

0. ( >45 

0 . 045 

0. 045 

0. 035 

II. ?o 

-o . ?< I 

l-.yu 

1 . < lO 

U . o52 

-0.015 

O. 045 

O. 043 

O. 047 

O. 03 7 

0 . 1 0 

—o . : o 

1 . 02 

1 .ol 

U. 040 

— O. Ol 2 

O. 043 

0. 04 3 

0. 043 

(>. (>?.7 

0. oo 

-u. JU 

1 . 07 

1 . 02 

( i . 07 7 

-0.008 

O. 048 

0.048 

0 . 050 

0.04? 

-o. to 

-o . “O 

1.06 

1 . 05 

O. 042 

-0. Ol 2 

0. 045 

0.045 

0. 045 

0. 042 

-O. 20 

- O- 20 

1 .09 

1 . 08 

0. 045 

-0 . f > 22 

V. 047 

O. 048 

0 . 04 7 

0. < >30 

-o.r.o 

-0 . 2‘ > 

1 . In 

1.10 

0.043 

-0 . 022 

0 . U50 

0 . 052 

0 . 053 

U. 040 

-O. 4(1 

-o. 20 

1.17 

1.17 

0 . 050 

-0. 023 

U .04 7 

0. 047 

O. 048 

o. 033 

-O. ‘in 

“0. Jo 

1.15 

1.15 

0.064 

-O. 023 

0.045 

0. 045 

0. 043 

0.038 

-0.55 

-0.'"O 

1.17 

1.16 

0. 0 72 

-0.028 

0.048 

0. 04 FI 

0. 04? 

0 . 037 

— U. 61 i 

-0- 20 

1.17 

1 . 16 

O. 070 

-0.022 

O. 045 

O. 047 

0. 045 

0 . 038 

-o. 65 

- O . /(") 

1 . IB 

1.17 

0. 074 

-0.0? 1 

O. 048 

0. 050 

0. 043 

1 1 . 038 

-0. 7u 

-0. ?o 

1 . 20 

1 . 20 

ft. 074 

-O .018 

0. 053 

0. 053 

0. 040 

< i . 04 5 

- 0 . 75 

~o. 20 

1.22 

1.21 

0.075 

-O.olo 

0 . 04B 

(>. 048 

0 . 038 

0.045 

-0 . 3o 

-0. ?o 

1 . 22 

1 . 22 

0. 077 

-0.0O2 

0. 048 

0. 048 

0. 040 

(i .04 7 
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y 

z 

Urms 

u 

V w 

U' rms 

u' 

v' 

w 1 

-O. GO 

-0. 30 

1.2;; 

1 . 22 

0.084 0.001 

0 . 055 

I). 057 

0.045 

0. (>42 

-O. 75 

— f> . 50 

i .?i 

1.21 

0.075 -0.013 

0 . 050 

0 . 050 

0.03.3 

0 . 040 

-O. 70 

-O. 50 

i.ro 

1.19 

0.077 -0.022 

0. 045 

0. 045 

0.04 2 

O. 043 

-0,63 

-0 - 30 

t . 18 

l . 18 

0.077 -0.028 

0. 047 

0. 047 

0. 045 

0. 043 

- 0 - 60 

-O. 30 

1.17 

1 . 16 

0.070 -0.028 

0 .048 

0. 048 

0.043 

O. 040 

-0. 35 

-0. 50 

i . jo 

1.16 

0. 072 -0. Olu 

0. 048 

0. 050 

O. 040 

<>. 042 

-0 . 5( > 

-0- 30 

1.15 

1.15 

0.070 -0.028 

0. 048 

0. 048 

O. 047 

(1. 042 

-0. 4o 

— 0 . 30 

1.12 

1. 12 

0.057 -0.010 

0. 047 

0. 047 

0.052 

0. 043 

— o. "Jo 

— O. 3( > 

1. 10 

1.10 

0.043 -0.002 

0. 028 

O. 028 

0.037 

0.022 

-0 . 20 

-0 . 3c i 

1 . 08 

1 . 07 

0.052 -0.020 

O. 050 

0. 050 

0. 057 

0. 047 

-O. 1 * ) 

-0. 5o 

1 . 05 

1 . 05 

0.047 — O .018 

0 . 050 

0. 048 

0 . 050 

0. 042 

O .(Ml 

-O. 3u 

t . <>3 

l.o3 

0.048 -0.013 

0. 047 

O. 047 

0 . 050 

0. 042 

o. lo 

- 0 - 3.0 

1 . U 1 

1 . Ol 

<1.042 -0.005 

0. 052 

0. 052 

0. 058 

O. 045 

O. 20 

-O. 3u 

0.99 

O. 99 

0.043 -0.012 

0. 047 

0.04 7 

0 . 050 

U. 043 

0 - 30 

-0 . "SO 

0. 96 

0. 96 

0.012 -0.017 

0. 043 

0.043 

0.048 

U. 038 

O. 40 

-0 . 30 

0.94 

0.94 

0 . 028 -0 . 030 

0. 043 

0.048 

0.048 

O. 042 

O. 50 

- o - 30 

U. 94 

0. 94 

O.02O -0.U27 

0. 047 

0.048 

0.048 

0. 035 

c > . 60 

-0-30 

0.91 

0.91 

0.015 — O .023 

0.048 

0.048 

0 . 052 

O. 035 

U. 70 

-c ) . 30 

0.88 

o. 81.1 

0.018 -0.022 

U. 063 

O. 063 

0.058 

0. 043. 

t » . BO 

-O. 30 

0.85 

0.85 

-0.003 -0.008 

0.060 

0 . 060 

0 . 062 

( i . 050 

o. 0 o 

- 0. 40 

It. 84 

0.84 

0.000 -0. 003 

O. 072 

U. 072 

U .060 

O. 045 

0 . 7 u 

—O. 40 

0. 89 

0 . 88 

0.007 — 0. U?7 

0 . 055 

0. 057 

0.057 

0. 047 

0 . 6< > 

-0. 4o 

O. 92 

0.91 

0.020 -0.032 

O. i )45 

0.047 

0. 04 7 

o.i >40 

o . 50 

-0- 4o 

o.?4 

0.91 
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0 . 03? 

0.027 

-0 . 05 

-0. 30 

1 . 24 

1 . 24 

O. 076 

-0. 016 

0.032 

0. 032 

0.044 

0 . 030 

-O. 05 

-U . 40 

1.71 

1.21 

0. 052 

-0. 035 

0.043 

0. 043 

0 .064 
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-0 . 85 

-0 . 50 

1 . 15 

1.15 

0.059 

-0.037 

0.078 

0. 079 

0 . 065 

0. 04 3 

-0 . 85 

-0 . 60 

1.15 

1.14 

O . OOo 

-O . 02 1 

0.059 

0 - 059 

0 . 059 

0.040 

-0. 85 

-O. 70 

1. 17 

1.17 

-0.022 

0.010 

0. 046 

0. 046 

0 . 059 

0. 038 

-0. 85 

-0 . 80 

1.10 

1.18 

-O. 043 

0.027 

0.037 

0 . 037 

(). 044 

0 . 032 

-0. 90 

0. 80 

1 . 08 

1.05 

-0.071 

-0. 244 

0. 103 

0. 104 

0. 091 

0. 083 

-0.90 

U. 70 

1 . 06 

1 .03 

0 . 005 

-0.226 

0.081 

0. 084 

0. 064 

0.076 

-0. 90 

0 . 60 

0. 98 

0.96 

0. 015 

-0. 202 

0. 008 

0 . 093 

0. 059 

ft. 083 

-(>. 90 

0 . 50 

U. 94 

0.92 

0. 03? 

-O. 160 

0.079 

0. 084 

0 . 059 

0.079 
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0. 40 

O. 82 

0.81 

O. 05? 

-O. 123 
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0. 098 

0. 066 

0 . 06 7 

-0. 9< t 
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0. 95 

0. 94 

0. 051 
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0. 145 

ft. 061 

0 . 066 
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1 . o3 

1 . 02 

0. 044 

-0 . <>57 
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0 . 06 1 
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0 . 1 U 

1 . 09 

1 .09 

O. 046 
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f>. 1 18 
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o. 052 
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1.13 

O. 044 
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0 . 057 
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0.96 
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0 . 06 1 
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0. 85 

0.82 
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0.071 
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0. 037 
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O. 083 
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0. 067 
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0.91 

-0. 002 
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0. 083 
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ft. 06? 
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-0 . 9( ) 

~U. 70 

1 . 00 

0.90 
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0. 088 
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0 .066 
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0 . 083 
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0 . 059 
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0.58 

-0.045 
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-0. 95 
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0.94 
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0. 073 

0 . 097 
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0 . 6( t 
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0. 159 
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0 .51) 
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0.88 

0 . 058 
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0. 067 
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0.77 
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0. 070 
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0 . 30 

(I. 79 

0.77 

0. 020 

-O. 139 
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0.85 

0. 023 
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O. 142 

0. 062 
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0.87 
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- 0 . 065 
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0.82 
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ft. 142 
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ft. 189 

0.259 
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0.93 
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O. 124 
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(t. 142 
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0 . 1 30 
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-U . 95 
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0.86 
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O. 094 

-0. 95 
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0 . 065 
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0.75 
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0. 1 19 
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0. 082 
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0. OOO 
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U. O7o 
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-0 . 40 

0 . 66 

0.64 
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0. 067 

0. 074 
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-0.50 
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0.81 
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0.86 
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0.114 

• O. 95 

- U. 70 

0.95 
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0.79 
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0.89 
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0. 140 
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0. 070 
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0.96 
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0.81 

0.71 
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O. 172 

0.062 

0.11 5 
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O. 72 
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0.070 
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O. 63 
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0. 12 7 
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0 . 50 
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O. 170 
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-0.013 
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0. 075 

0. 075 
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-o. to 

O. 44 
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-0.002 

0. 027 

0.21 0 

0.21? 

0.067 

0. 075 

-0.96 
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0. 59 

0.50 

0 . OOO 

0.058 

0. 133 
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0 . 058 

0.U7B 

-0. 96 

-0 . 30 

O. 63 

0. 61 

0 . 000 

0.117 

0.112 

0. 1 15 

0 . 055 

0. 073 

-0.96 

-0. 40 

0 . 58 

0 . 56 

0.007 

0. 123 

O. 105 

O. Hu 

0. 068 

0.070 

-0. 96 

-0. 50 

0.71 

0.67 

0. 027 

O. 210 

0. 142 

0.144 

0.053 

0 . i >93 

-O. 96 

-O. 6< > 

0.86 

0.81 

0. 002 

0. 254 

0. 144 

0. 157 

0 . 058 

0.O07 

-0. 96 

— O. 70 

0.97 

0.81 

-O. 012 

O. 290 

O. 150 

0. 175 

0. 068 

0 . 005 

- U. 96 

- 0 . BO 

0.90 

0.95 

-0. 025 
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0. 1 32 

u. 139 

O. 062 
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V 
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u' 

v' 
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0. O0 

1 . 00 

0.99 

O. 065 

-O. 020 

0. 063 

0.058 

0. 122 

O. 048 

O. 70 

0. 00 

1 . 03 

1.02 

0.058 

-O. 005 

0.055 

0.053 

0. 082 

0.042 

0 . 60 

0.00 

1 . 05 

1.04 

0. 072 

0.000 

0.050 

0. 048 

0.063 

0. 033 

0.50 

0.00 

1 . 05 

1.05 

O. 080 

0.000 

0.042 

O. 040 

0. 060 

0 . 030 

0. 40 

O. 00 

1.06 

1 . 06 

0.090 

0.007 

0. 035 

0.035 

0.037 

0.023 

0.30 

O. 00 

1 . 08 

1.07 

0. lOO 

0.007 

0. 040 

0.040 

0.028 

0.025 

O. 20 

0.00 

1.09 

1.08 

0.112 

0.010 

0.038 

0.038 

0.032 

o. 027 

O. 10 

0 . 00 

1 . lO 

1.09 

O. 122 

0.008 

0 . 037 

0.038 

0.030 

0.018 

0. 00 

0 . 00 

1.11 

1.10 

0-1 33 

0.008 

0.038 

0.037 

0. 038 

0 . 023 

-0. 10 

0 . 00 

1.13 

1 . 12 

O. 153 

0.010 

0.040 

0. 040 

0.038 

0.023 

-V. 20 

O. Ou 

1. 14 

1.13 

O. 167 

0.008 

0. 040 

0. 040 

0. 032 

0. 027 

-0. 30 

o . oo 

1. 15 

1. 13 

O. 1B8 

O. 010 

0. 040 

0.040 

0.035 

U. 035 

-O. 40 

0 . oo 

1 . 16 

1.15 

O. 205 

O. 013 

O. 040 

0.040 

0* 033 

0. 023 

-0.50 

i » « 00 

1. 18 

1.15 

0. 235 

0.013 

0. 043 

0.043 

0. 042 

U - 030 

-0. 55 

0. 00 

1 . 16 

1. 13 

0. 238 

0.020 

0.070 

0. 072 

0.050 

0.053 

-0.60 

0 . 00 

1.05 

1.02 

0. 223 

0.017 

0. 140 

0. 1 43 

0.063 

0.082 

-0.65 

0 . < >0 

0. 86 

0. 83 

0. 197 

-0.005 

0.138 

0. 142 

0.078 

0.090 

-0.70 

0. 00 

0. 73 

0. 70 

O. 158 

-0.013 

0.113 

0. 115 

0. 072 

o. 092 

-0.75 

0. 00 

0.66 

0. 64 

0. 130 

-0.053 

0. 093 

0.093 

0. 060 

0. 087 

-0 . HO 

0.00 

0.64 

0.62 

0.090 

-0.077 

0.093 

0.092 

0.058 

0. 09? 

-0. 80 

-0. 1 0 

0.60 

0.59 

0.077 

0.027 

0.093 

0 . 093 

0.055 

0.077 

-O. 75 

-0. 10 

O. 65 

0.63 

0. 118 

0.018 

0. 1 OO 

0. 10? 

0. 062 

0. n85 

-0. 70 

-o. lo 

O. 75 

0. 73 

O. 153 

0- 000 

0. 123 

0. 125 

0.070 

0 . OR3 

-0.65 

-0. lo 

0. 86 

0. 83 

0. 183 

-0 .018 

o. 132 

0. 137 

0.O73 

0 . 092 

-O. 60 

-0. 10 

1.05 

1 . 02 

O. 220 

-0- 020 

0. 132 

0. 137 

0.073 

0. 075 

- 0 - 55 

-o . 1 0 

t. 15 

1. 12 

0- 230 

-0.017 

O.067 

0 . 0 70 

U . 057 

0.051 

-0. 50 

-0 . 1 o 

1.17 

1.15 

U. 230 

-0.008 

0 . 04 7 

0. 04 7 

0 . 038 

0. 02 3 

- (». 40 

“0. to 

1. 16 

1.14 

0. 203 

-O . 003 

0. 040 

0. 040 

0.035 

O. 022 

— 0 . 3n 

-0. 10 

1. 15 

1.13 

O. 180 

-0.005 

0. 037 

0. 038 

0.035 

0. u20 

-<>. 20 

-0. 10 

1. 13 

1.12 

0. 162 

0.000 

O. 040 

0. 040 

0. 033 

0. 020 

-0. 1 0 

-0. 10 

1.12 

1.11 

0. 145 

-0.002 

0. 03B 

0 . 038 

0 . 030 

0 . 022 

0. 00 

-0 . 1 o 

1 . 1 1 

1. 10 

0. 137 

0.000 

0 . 038 

0 . 038 

0 . 033 

U . 020 

0. 10 

-0. 10 

1.11 

1 . lo 

0. 122 

0.002 

0. 040 

0. 040 

0. 030 

0.023 

0. 20 

-0. 10 

1 . 08 

1.07 

0. 112 

-0 . 002 

0. 042 

0. 042 

0.032 

0. 027 

0 . 30 

-0. lo 

1.07 

1 . 06 

0. 102 

-O. 005 

0.042 

O. 042 

0.035 

0.025 

0 . 40 

- 0 . to 

1.07 

1 . 07 

0.090 

0 . < >00 

0. 042 

0.043 

0 . 035 

U. 025 

0 . 50 

-0 . 1 0 

1 . 05 

1 . 05 

0.077 

-0.002 

0. 042 

0. 042 

0.035 

0.025 

0 . 60 

-0. lo 

1 . 03 

1-03 

0.062 

-O . 007 

0.047 

0.047 

0.045 

0 . 038 

0. 70 

-0. 10 

1 . ul 

1 . ul 

0 . 055 

-0.012 

O. 053 

0.053 

0. 057 

0. 045 

<>.8o 

-0. lo 

0.98 

0. 97 

0. 053 

-0.013 

0.062 

0 . 065 

0 . 1 02 

0 . 052 

0. 00 

-0. 20 

1 . 00 

0.99 

0. 055 

-0.015 

0.070 

0 . 068 

0. 105 

0. 057 

O. 70 

-0. 20 

1 .01 

1 . 00 

0. 057 

-0.012 

n. 053 

0.053 

0.092 

0.040 

0. 60 

-O. 20 

1 . 04 

1 . 04 

0.075 

-0.013 

0.048 

0.048 

0. 050 

0. 03 7 

0 . 5*> 

-0. 20 

1 . 05 

1.05 

0.077 

-0.007 

0 . 043 

0.043 

0.043 

0 . 030 

0. 40 

-0. 20 

1 . 06 

1 . 05 

0. 088 

-O . 005 

0. 040 

0. 040 

0 . 035 

0.023 

0 . 30 

-0. 20 

1 . 08 

1 . 07 

0. 102 

-0.018 

0. 042 

0. 042 

0.033 

U . 0?R 
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1 . to 

1 . 09 

0. 1 17 

-0.017 

0. 042 

U. 042 

0. 043 

0 - 030 
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-0. 20 

1.11 

1.10 

0. 127 

-0.015 

O. 038 

0.038 

0. 045 

0. 02 7 
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-0. 20 

1.12 

1.11 

0. 133 
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O. 038 

0 . 038 

0 - 050 
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-o . 20 

1.12 

1.11 
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-0.015 

o. 040 
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0.045 
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-o. 20 

-0. 20 

1 . 15 

1 . 13 

O. 160 

-0.018 

0. 042 

0.042 

0.043 

0.025 

— 0. 30 

-0. 20 

1.15 

1.14 

0. 177 

-0.025 

O. U42 

0. 042 

0 . 035 

i ) . ul'3 

-0. 40 

-0 . 20 

1.17 

1.15 

0. 195 

-0.037 

O. 043 

0. 045 

0 . 033 

0 . 025 

-0- 50 

-0. VO 

1 . 18 

1.16 

0. 213 

-0.050 

0.048 

0.047 

0. 048 

O. 02 7 

-0.55 

-0. 20 

1 . 17 

1.15 

0. 223 

-0.058 

0.062 

0.063 

0. 048 

0 . 030 

-0 . 60 

-0 . ?t ) 

1 .04 

1-02 

0.210 

-0.048 

0. 132 

o. 135 

0 . 068 

0 . 063 

-0.65 

~l>. 20 

0. 88 

0.85 

0. 190 

-0.043 

0. 127 

0. 130 

0.072 

0. 080 

-0. 70 

“ 0 . 20 

o.7a 

0.75 

0. 173 

-0.U03 

O. 098 

0. 100 

0. 072 

O. 0 73 

-0. 75 

-0 . 2u 

O. 72 

0. 69 

O. 153 

0.038 

0. 092 

O. 092 

0 . 063 

0. 077 

-O.8o 
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0.64 

0.123 

0.073 

i >. ( >98 
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0 . 068 
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0.87 
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-0. 30 

0.87 

0. 81 
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-O. 033 
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0 . 068 
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-O . 30 

0.91 

0.89 
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\PPENDIX » 

High Reynolds Number Tabulated Data 


Table D-l Station 2, X(13) = 2.0 Fu trance 
Table D-2 Station 4, 0 - 0° 

Table D-d Station G, 0 = 30° 

Table D-4 Station 8, 0 = GO 0 
Table D-r* Station 10 0 = 00° 

3 able D-C Station 11, X(H) = 1.0 Exit 


The funnel coordinate positions have been nc n-dimcusionaU’Acc on dnet half- 
3 >li //.. A.b ■vekeHies have been noM-dimcnsiooaUy.od on bulk ve’ocity (20 m/s). 
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- C>. 0u8 

o. 029 

0.029 

n. 021 

u.iilB 

0.75 

0.00 

u . 95 

0.95 

- o . 1 105 

-0 . 005 

Ci . u30 

O. 030 

0. 020 

0.018 

0.70 

0 . 00 

o. 95 

0.95 

- 0. 1108 

- 0.000 

0.031 

O. 031 

U .019 

0.017 

0. 65 

o . oo 

O. 95 

0.95 

-0.01 1 

-0. 008 

0.029 

0.029 

0. 019 

o. oi a 

V . 60 

o . ao 

0 . 95 

U . 95 

0. 0 1 5 

-0. ol 1 

0.029 

0. 029 

0 . 1 120 

0.017 

0. 55 

0. Oo 

0. 95 

0. 75 

-o . 0 i a 

-0.010 

0. o 2 9 

0.029 

U . 020 

0 . 0 1 6 

0 . 50 

0 . Ro 

0 . 9 4 

0.94 

-0.1114 

-0.012 

0.027 

0. 02! ! 

Cl . U ] 7 

ii. 017 

0. 45 

O. 00 

O. 75 

0.94 

Cl. 014 

-0. Lit 1 

0.027 

Cl. 027 

O. 01 B 

0 . 0 1 7 

O. 40 

o . uo 

0 . 95 

O. 95 

-0.013 

-0. 007 

0.020 

0 .028 

0. 020 

U. 017 

0. 30 

o. ao 

0. 95 

0. 75 

-0. 01 1 

-0. OuB 

0.027 

0.027 

0.019 

Cl. ol7 

0. 20 

o . ao 

0. 74 

0.94 

-0. 007 

-0.010 

0.0 33 

0 . 033 

0. 018 

0 .019 

0. 10 

0 . 00 

0 . 70 

0. 90 

-O. 003 

0. 003 

0.029 

0. 029 

O. 023 

0. 020 

0 . 00 

0.80 

0.91 

0.91 

-O. Ci 12 

Cl . 0u2 

0.024 

0. 024 

O. 02 3 

0.019 

-0. 05 

o. ao 

0. 93 

O. 93 

-0.015 

-O. 007 

O. 020 

0. 020 

0. 022 

0. 016 

- 0 . lo 

0 . 80 

0.94 

0.94 

-0. 014 

-Ci. 008 

0.022 

0 . 022 

0. 023 

0 .016 

-0. 15 

0 . BO 

o . 94 

0.94 

-0 .011 

-0.009 

0.025 

0 . 025 

0.021 

0 .017 

-0. 20 

o. ao 

0. 94 

0. 94 

-0. 006 

-Cl. Clio 

0. 028 

0.028 

0.017 

0 .017 

-0. 25 

0 . 00 

0. 95 

0 . 95 

-0.007 

-Ci. 006 

0.030 

0 . Ci 3 0 

0 . 020 

0.018 

-0. 30 

0 . Do 

0.95 

0 . 95 

-0. U03 

-0. 002 

0. 032 

Ci. 032 

Ci. 019 

0.018 

-0 . 35 

0 . 80 

0.97 

0. 97 

-0. 0u5 

-O. 003 

0. 035 

0. 035 

0. 020 

0 .018 

0 

1 

U.BO 

0. 79 

0.99 

-U. 006 

0.1102 

0 . 034 

0. 034 

Ci . 020 

0. 0113 

-0. 45 

0 . So 

1 . 00 

1 . Oo 

-0. Li Cl 4 

0. 005 

0 . 03 1 

0 . 03 1 

0.021 

0. 1.115 

-0. 50 

0. BO 

t.UO 

t . Oo 

-0. 007 

O. 005 

0. U29 

0.029 

0. 022 

ci. 010 

-0.55 

o. ao 

0.99 

0.99 

-0.004 

0.002 

0 . 03 2 

0. 032 

0.017 

O. 01 0 

-0. 60 

0.80 

1 . 00 

1 . Ou 

-0. 004 

0. 004 

0.034 

0. 034 

0. 019 

0.O16 

—0. 65 

o. SO 

1 . oO 

1 . 00 

-0. 005 

O. 004 

0. 028 

Ci . 023 

0 . 020 

0 .017 

-O. 70 

o. ao 

1 . 00 

1 . in') 

- 0. 002 

O. 004 

0.029 

0.029 

0.019 

C> . i i 1 9 

-0 . 75 

0. 00 

1 . 00 

1 . oo 

-0.001 

0. 005 

O. 029 

0. 029 

0 .019 

Ll. 010 

-o . ao 

O . Bo 

1 . 00 

1 . 00 

0.004 

0 . 003 

0 . 026 

0.026 

0.018 

u.015 
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_y 

Z 

Urms 

U V 

j 

1 

U* rms 

u' 

vi 

w' 

o.ao 

0* 00 

1.03 

1.03 -0.002 

-0.007 

0.030 

0.031 

0.023 

0.021 

O. 70 

0. ou 

1.04 

1.04 -0.019 

-0. 004 

0 . 033 

0. 033 

0.024 

0.020 

0. 60 

0 . 00 

1.04 

1 . 04 -0. 028 

-0.007 

0.033 

0.034 

0. 018 

0. 020 

0 . 50 

U . 00 

1.05 

1.05 -0.033 

-0 . 002 

0.031 

0.031 

0.019 

0.021 

0. 40 

0. 00 

1 . 05 

1.05 -0.034 

-0. 009 

0.038 

0. 039 

0.019 

0.022 

0.50 

0.00 

1.05 

1.05 -0.027 

-0.014 

0. 036 

0.037 

0.021 

0.023 

0.20 

0. 00 

1 . 06 

1.06 -0.033 

-0. 007 

0.041 

0. 041 

0. 020 

0.021 

0. 10 

0 . ijQ 

1.06 

1.06 -0.031 

-0.01 1 

0.041 

0.041 

0.020 

0.021 

o . no 

0. 00 

1.07 

1.07 -0.030 

-0.009 

0. 040 

0. 041 

0.019 

0.021 

-0. 10 

0 . u 0 

1.07 

1.07 -0.031 

-0.008 

0.042 

0.042 

0.019 

0.021 

-0. 20 

0. 00 

1 . 08 

1 . 08 -O. 029 

-0.014 

0. 040 

0. 040 

0. 020 

0. 020 

-0.50 

0. 00 

1 . 09 

1.09 -0.027 

-0.010 

O. 042 

0. 042 

0. u20 

0.019 

-0. 40 

0 . 00 

1 . 09 

1.09 -0.026 

-0 .013 

0. 043 

0. 043 

0 .018 

0 . 02 1 

-U.50 

0 . oo 

1.10 

1.10 -0.029 

-0 . 009 

O. 044 

0. 045 

0 .017 

o. c>2l 

-0. 55 

0. 00 

1.11 

1.11 -0. 027 

-0.01 1 

0. 043 

0. 043 

0.019 

0 . 020 

-U. 60 

0. 00 

l. 1 1 

1.11 -0.027 

-0.01 1 

O. 042 

0. 042 

0.1(17 

0. U 2 1 

— 0 . 65 

0 . 00 

1 . 12 

1.12 -0.025 

-0.01 1 

0.041 

0 . 04 1 

O. 017 

0.023 

-0. 70 

0 . uO 

1.13 

1.13 -0.021 

-0.010 

0. 041 

O. 041 

0 .013 

0.022 

- 0. 75 

0. 00 

1 . 13 

1.13 -0.019 

-0 . 008 

0 . 04 1 

0 . u4 1 

0.017 

(1.021 

— U. 80 

0. 00 

1.14 

1.14 -(). 01 6 

-0.008 

0 . 04 2 

0. 047 

0 .016 

0.021 

-o. no 

“0. 1 0 

1.14 

1.14 -0.014 

-0.010 

0. 041 

O. 04 1 

0. u20 

0. 023 

-0.75 

-0.1 0 

1.13 

1.13 -0.017 

-0.011 

0. 040 

0. 04 0 

0.019 

0. 021 

-O. 70 

- 0 . 1 o 

1 . 13 

1.13 -0. 022 

-0.010 

0. 040 

0. 040 

0 .019 

0. 02 1 

-O. 65 

“ 0 - to 

1 . 12 

1.12 -O. 023 

- 0. 008 

0. 040 

0. 040 

0.019 

0. 022 

-<>. 60 

-0. 10 

1 . 12 

1.12 -O. 029 

-0 . 009 

0.04 1 

0. 041 

0 . 030 

U. 022 

-0.55 

“O. to 

1.11 

1.11 -0.030 

-0.008 

U. 044 

0.044 

0 .018 

0 .021 

-0. 50 

-0. 10 

1.11 

1.11 -0.030 

-0. 009 

0. 046 

0.046 

0.071 

0 . i.i" 1 

• o. 40 

~o. io 

1 . iu 

1.10 -0. 02o 

-O.olO 

0. 048 

0. U48 

0 .019 

0 . 0 20 

-0. 30 

-0. 10 

1 . 10 

1.10 -O. 026 

-0.012 

0. 044 

0. 044 

0 . 02 0 

0. 022 

-- 0 . 2 O 

-o.l 0 

1 . 08 

1.08 -0.027 

-0.014 

0. 04 3 

0. 043 

0 . 020 

0 . 073 

-0. 10 

-0. lo 

1 . 07 

1.07 -0.027 

-0.012 

O. 046 

0. 046 

0 . u22 

0. 021 

0. 00 

0 . 1 0 

1 . 07 

1.07 -0.029 

-0. 010 

0. 079 

0.039 

0. 020 

0. 021 

0.10 

-0 . 1 0 

1 . 06 

1.06 -0.032 

-0.013 

0.038 

0. 033 

0.030 

0 . u2 1 

O. 2(1 

-0. 10 

1 . 05 

1 . 05 -O'. 037 

-0.O0B 

( i . 034 

0 . 0 35 

0 .019 

0. 071 

0. 30 

-0. 10 

1 . 05 

l.uS -0.031 

-0.010 

0 . 035 

0. 035 

0 . 020 

0 . 022 

0. 40 

"0. 10 

1 . U5 

1.05 -0.036 

-0. 004 

U. 035 

u. 035 

0.018 

0 . 02 1 

0 . 50 

-o. 10 

1 . 05 

1 . 05 -0 . 034 

-0.006 

0 . 037 

0 . 037 

0.019 

0. 022 

0. 6o 

0. 10 

1 . 04 

1.04 -0.032 

• 0.003 

O. 033 

0 . < »33 

0. 070 

0 . 0 1 9 

0. 70 

-0. 10 

1 . 04 

1.04 -(J.023 

-O. U05 

0 . ('<3 2 

0. 032 

O. 025 

0. 02? 

0 . 80 

-O. 1<> 

1 . 03 

1.03 -0.010 

• 0. 004 

0.0 30 

u . O 30 

(>. 077 

0.((?1 

o . no 

- 0. Zo 

1 . 03 

1 . 0.3 -0.014 

-0. 001 

0. 032 

0. 032 

o. 026 

((. 02? 

0. 70 

■ On 20 

1.04 

1.U4 -0.017 

-u.oul 

0. 030 

<i. o30 

0. 024 

0.H21 

0 . 60 

-0 . 2< > 

1 . 04 

1 . u4 -0.029 

-0. (105 

u . 034 

0. 034 

0 . 020 

0. ((22 

0 . 50 

-0. Z'u 

1 . o5 

1.05 -0.033 

-O . 002 

U. 032 

U. 02.2 

0. ((20 

t'(. 022 

0. 40 

— 0 . zo 

1 . 05 

1.05 -0.035 

-0. 007 

0.037 

0. 037 

0. 020 

0 .023 

0 . 30 

- o . :*o 

1 . 05 

1 . 05 -0.035 

- o. 007 

0 . 038 

0.U38 

o . or i 

0 . 02 1 

O. 20 

— 0 . zo 

1 . 06 

1 . (i6 -O. 036 

-0. 007 

O. 039 

0. 039 

0. 021 

0. 021 

0 . 1 0 

- 0 . 2o 

l . 06 

1.06 -0.U32 

-O. 012 

0. 037 

(>. 038 

0. 021 

0. 073 

0. 00 

-0. 20 

1 . 06 

1 . u6 -0.031 

-0. 01 1 

U. 042 

0. 042 

0 . 020 

0. o21 

•0. 10 

- 0. 20 

l.(»7 

1.07 -0.034 

-O . 008 

0 .037 

0. (>39 

0. 021 

O. 022 

-0. 70 

-0 . zo 

1 . 08 

1.08 -0.032 

-0. 01 4 

O. 043 

0. 048 

0. 019 

0 . 02 1 

o 

r- 

0 

1 

- u. 20 

1 . 09 

1.09 -0.033 

-0.013 

0. 045 

0. 046 

0. 02( ( 

0. 022 

-0. 40 

-0. 20 

1 . 10 

1.09 — u. 034 

-O . 0 1 2 

0.04 1 

0.041 

0 . 022 

0. 020 

- 0 . £»0 

0 . zo 

1.11 

1.11 -0. 033 

-O. 015 

0. 044 

0.044 

0. 019 

0. 020 

-n. 55 

20 

1.11 

1.11 -0.030 

-0 .012 

0. 045 

0. 045 

0. 018 

0. 022 

-0 - fc>u 

~0- 20 

1.11 

1.11 -0.023 

-0.013 

0. 041 

0.04 1 

0 .019 

0. 023 

* 0. t»5 

-0- 20 

1.11 

1.11 -0.025 

-O. 012 

O. 040 

0. 040 

0. 020 

0. 022 

-0. Hr 

-II- Zo 

1.12 

1.17 -0.027 

-O. 014 

O. 038 

0 . 038 

0.019 

0. 022 

- 0. 75 

-o. ZO 

1.13 

1.13 -0.023 

-O. 009 

0. 038 

0. 038 

0. 019 

0. 022 

— U . HU 

-0 . 20 

1.14 

1.14 -0.021 

-O . o 1 0 

0 . U4 1 

0.041 

0.019 

0. 021 
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_ _y 

z 

Urms 

JL_ 

V 

w 

U' rms 

u' 

v' 

w ' 

-O. 80 

~o. 30 

1.14 

1. 1.3 

-0. 023 

-0 . 008 

0 . 039 

0. 039 

0.017 

0. 022 

-0 . 75 

-0. 30 

1 . 13 

1 . 13 

-O. 075 

-0.012 

0. 040 

0. 04 1 

Ci. 019 

0. 022 

-0. 70 

--Ci. 3c i 

1.13 

1.13 

-O. 025 

-0.012 

0. 042 

0 . 04 2 

0.018 

0. 022 

-O. 65 

-0. 3<"> 

1.13 

1.13 

-O. 023 

-0. 01 1 

0. 1544 

0. 044 

0. 021 

0. 024 

-0. 60 

-0 . 3 Ci 

1 . 12 

1.11 

-0.027 

-0. 1512 

0.045 

0 . 045 

0. 020 

0. 020 

-o . 55 

-0 . 30 

1.10 

1 . 1 Ci 

-Ci. 028 

-0.013 

0.042 

O. 042 

0. 021 

0 . 022 

-0. 50 

-0 . 30 

1 . lo 

1.10 

-0.033 

-0.012 

0 . 043 

0. 048 

0.019 

0. 021 

-0. 40 

-0. 30 

1. 10 

1.10 

“0.033 

-Ci. CilO 

Cl. 046 

Cl. 046 

0. 021 

n. 020 

-0. 30 

-0. 30 

1 . 09 

1.08 

“0. 030 

-0. 008 

0. 048 

0. 048 

0. 022 

0. 021 

-0.20 

- 0 . 30 

1.08 

1 . 08 

-0. 029 

-0. 017 

0. 046 

0.047 

Ci. 1571 

0.1 >22 

-0. 10 

— o. 30 

1 . 07 

1 . 07 

-0. 028 

-O. c iO? 

0. c544 

0. 044 

0.019 

0. 032 

o . oo 

-0. 30 

1 . 06 

1 . 06 

-0.025 

-0.010 

0.039 

O. 039 

0 . 02 1 

0. 022 

0. 10 

-O . 30 

1 . 06 

1 . 06 

~ U. 027 

-0.01 1 

0.040 

0. 040 

0.019 

0. 022 

0 . 20 

-0. 30 

1 . 05 

1.05 

- 0 . 030 

-G.OlU 

0 . 039 

0.040 

0. 019 

0. 1522 

0. 30 

-0. 30 

1 . 05 

1 . 05 

-0. 030 

-O. 008 

0. 039 

0. 039 

0.021 

0. 022 

0. 40 

--O.30 

1 . 04 

1 . 04 

-0.028 

-0.006 

0 . 034 

0. 034 

0.020 

0.022 

0. 50 

-O. 30 

1.04 

1.04 

-0.031 

0 . 002 

0.035 

0.035 

0. 018 

0. 021 

0.60 

-Ci. 30 

1 . 04 

1 . 03 

-0.025 

-0 . 006 

0. 034 

0. 034 

0.021 

0. 02 Cl 

0. 70 

-Ci . 30 

1 . 03 

1 .03 

-0.016 

-0. 008 

0.034 

0.034 

0.026 

0. 022 

0. 80 

-0. 30 

1 . 03 

i.o’; 

-0.016 

-O. 005 

0. 029 

0.029 

O. 026 

0. 020 

0 . 80 

-0. 40 

1 . 02 

1 . 07 

-Ci. 017 

—0.007 

0. 029 

0. 029 

0.024 

0. 1522 

U. 70 

--o. 40 

1.03 

1.03 

-0. 023 

-0 . O 1 o 

0.031 

0.032 

0. 029 

0 . 02 1 

ci. 60 

-Cl. 4c 1 

1 . c.i 7 

1 . 03 

-0. 022 

-0. 01 1 

0. 035 

0.035 

0.027 

0 . 07 1 

0 . 50 

--Cl ■ 40 

1 . 04 

1.04 

-1.1. 028 

-Cl. 012 

0.037 

0.037 

0. 072 

0.U22 

0. 40 

-0. 40 

1 . 05 

1 . 03 

-0. 027 

-O. 013 

0. 038 

0. 038 

0.019 

O. 021 

O. 30 

-0 . 4c i 

1 .05 

1 . 1)5 

-Ci. 028 

-0.013 

O. 040 

0. 040 

0 . 1520 

0.071 

O. 20 

— c l. 4 Cl 

1 . 06 

1 . 06 

-U. 030 

-0.01 1 

0.038 

0.038 

Cl. 021 

U. 019 

0. Id 

0. 40 

1 . oo 

1 . 06 

-0.029 

-0 . o 1 0 

Cl. 040 

0. 040 

0. 023 

O. 022 

( ) . i >0 

-Cl. 40 

1 . 07 

1.07 

-<i. c 130 

- Cl. 013 

0 . 04 1 

0. 041 

0.021 

0. 072 

-o. ill 

Cl. 4u 

1 . 07 

1 . 07 

- 0 . 032 

-0. 010 

0. ' i43 

0. 043 

Ci. 023 

0. 021 

— ci. 20 

-O. 40 

1 .08 

1.00 

”0.035 

-0.012 

0.049 

o . 050 

0. 022 

0 . Ci 2 2 

-0 . 30 

-Cl. 40 

1.09 

1.09 

- c.i. 035 

-0.011 

0. 048 

0. 048 

0.021 

0.022 

-0. 40 

- 0. 4o 

1. 10 

1. id 

-0 . u33 

-O. 012 

0.044 

0. 044 

0. 020 

O. 021 

0 . 50 

- O. 40 

1 . 10 

1. 10 

-0. 028 

-0.014 

0.042 

0. 042 

0. C)?l 

O.C>7ii 

-0. 55 

-0.40 

1. 12 

1 . 12 

-0.026 

-O . 0 1 2 

0.047 

0.04/ 

0. 021 

0. 020 

-0. 60 

-0. 40 

1.11 

1.11 

-0.025 

-11.013 

0. 043 

u. 04 Z 

Cl. 020 

Cl. 022' 

—0. 65 

-0. 4<J 

1.12 

1 . 12 

-O. Ci 25 

-O. 014 

0.047 

0.047 

0. 020 

0 . 022 

-0. 70 

-- Cl. 40 

1. 12 

1 . 17 

-0. 02o 

-0. 009 

0. 040 

0. U40 

Ci. 02 1 

0 . 02 1 

--0. 75 

-0. 40 

1.13 

1 . 12 

-0. 023 

-0.01 1 

O. 040 

0. 040 

0. 018 

0.022 

“0.80 

-0 . 4c.) 

1. 13 

l. 13 

— <‘>. 020 

-0 . Ou9 

0 . 1543 

0. 043 

0. 020 

0 . 022 

-0. 80 

-0 . 50 

1.13 

1.13 

-0. 021 

-0 . 007 

<5. 041 

O. 041 

O. 020 

0. 023 

-0. 75 

-0.50 

1.12 

1. 12 

-0. 020 

-Cl. 009 

0. 041 

0.041 

0.021 

0.022 

— O. 70 

-0. 50 

1 . 13 

1. 12 

-13. 023 

-0- 013 

0.042 

0. 042 

0. 020 

0. 072 

-0.65 

- 0 . 50 

1.12 

1.12 

-0.024 

-0.010 

0.045 

Cl. 045 

Cl. 018 

0.022 

-0. 60 

- 0 . 50 

1.11 

1.11 

-0.027 

-0.010 

0.043 

0. 043 

0.019 

0. 023 

-0. 55 

-0. 50 

1.11 

1.11 

-0.02/ 

-0 .010 

0.044 

0 .044 

0.021 

o . 022 

-0. 50 

-0. 50 

1.10 

1 . to 

-0.027 

-0. 01 1 

0. 045 

0. 045 

0.O23 

0 . 022 

-0. 40 

-0. 50 

1 . 09 

1 . 09 

-0.029 

-CI. 010 

0. 046 

Cl. 046 

O. 027 

0 . Ci 2 2 

— 0 . 30 
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